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We describe a measurement of the W boson mass mw using 200 pb _1 of ^=1.96 TeV pp 
collision data taken with the CDF II detector. With a sample of 63,964 W — + ev candidates and 
51,128 W -> [xv candidates, we measure m w = [80.413 ± 0.034(stat) ± 0.034(sys) = 80.413 ± 0.048] 
GeV/c 2 . This is the single most precise mw measurement to date. When combined with other 
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measured electroweak parameters, this result further constrains the properties of new unobserved 
particles coupling to W and Z bosons. 

PACS numbers: 13.38. Be, 14.70.Fm, 13.85.Qk, 12.15. Ji 



I. INTRODUCTION 

The discovery of the W and Z bosons in 1983 [l[ 
confirmed a central prediction of the unified model 
of electromagnetic and weak interactions 0. Ini- 
tial W and Z boson mass measurements verified the 
tree-level predictions of the theory, with subsequent 
measurements probing the predicted 0(3 GeV/c 2 ) 
[3, 0] radiative corrections to the masses. The current 
knowledge of these masses and other electroweak pa- 
rameters constrains additional radiative corrections 
from unobserved particles such as the Higgs boson 
or supersymmetric particles. These constraints are 
however limited by the precision of the W boson mass 
mw, making improved measurements of mw a high 
priority in probing the masses and electroweak cou- 
plings of new hypothetical particles. We describe in 
this article the single most precise mw measurement 
5] to date. 

The W boson mass can be written in terms of 
other precisely measured parameters in the "on-shell" 
scheme as 



m w 



If 



■kolem 



c V2G F (1 - m 2 w /m 2 z ){\ - Ar) 



(1) 



where oiem is the electromagnetic coupling at the 
rcnormalization energy scale Q = mzc 2 , Gp is the 
Fermi weak coupling extracted from the muon life- 
time, mz is the Z boson mass, and Ar includes all 
radiative corrections. Fermionic loop corrections in- 
crease the W boson mass by terms proportional to 



\n(mz /rrif) for mj <C mz [4], while the loop con- 
taining top and bottom quarks (Fig. [1]) increases mw 
according to @ : 



Ar tb 



'SGpm^r 



ti 3 8V2vr 2 (r 



*w) 
2m 2 m 2 



ln(m 2 /m 2 ) 



(2) 



where the second and third terms can be neglected 
since m t ^> Higgs loops (Fig. [2]) decrease mw 
with a contribution proportional to the logarithm of 
the Higgs mass (m#). Contributions from possible 
supersymmetric particles are dominated by squark 
loops (Fig. [3]) and tend to increase mw- Gener- 
ally, the lighter the squark masses and the larger the 
squark weak doublet mass splitting, the larger the 
contribution to mw- The total radiative correction 
from supersymmetric particles can be as large as sev- 
eral hundred MeV/c 2 0]. 

Table U @ shows the change in mw for +lcr 
changes in the measured standard model input pa- 
rameters and the effect of doubling mn from 100 
GeV/c 2 to 200 GeV/c 2 . In addition to the listed 
parameters, a variation of ±1.7 MeV/c 2 on the pre- 
dicted mw arises from two-loop sensitivity to a s , e.g. 
via gluon exchange in the quark loop in Fig. [TJ The- 
oretical corrections beyond second order, which have 
yet to be calculated, are estimated to affect the mw 
prediction by ±4 MeV/c 2 [|. 
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FIG. 1: The one-loop contribution to the W boson mass 
from top and bottom quarks. 



The uncertainties on the mw prediction can be 
compared to the 29 MeV/c 2 uncertainty on the world 
average from direct mw measurements (Table HI)) . 
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FIG. 2: Higgs one-loop contributions to the W boson 
mass. 




FIG. 3: One-loop squark contributions to the W boson 
mass. 



which include results from four experiments, ALEPH 
[3, DELPHI [H, L3 and OPAL (3, study- 
ing y/s — 161 — 209 GeV e + e~ collisions at the 
Large Electron Positron collider (LEP), and from two 
experiments, CDF and D0 [13, GJ], studying 
yfs = 1.8 TeV pp collisions in Run I of the Fermilab 
Tevatron. The current experimental mw uncertainty 
is a factor of two larger than the uncertainty from 
radiative corrections, excluding the Higgs contribu- 
tion (Table Q| . The Higgs mass constraint extracted 
from the W boson mass is thus limited by the direct 
m\y measurement. The precise mw measurement de- 
scribed in this article has a significant impact on the 
world average mw ■ 



II. OVERVIEW 

A measurement of mw at a pp collider [2l| is com- 
plementary to that at an e + e~ collider. Individual u 
(d) quarks inside the proton can interact with d (u) 
quarks inside the anti- proton (or vice versa), allow- 
ing single W + (W~) boson production, which is not 
possible at an e + e~ collider. In addition, pp collid- 
ers have higher center of mass energies and W boson 
production cross sections. This provides high statis- 
tics for the leptonic decays of the W boson, which 
are studied exclusively because of the overwhelming 
hadronic-jet background in the quark decay channels. 



Parameter Shift 


mw Shift 




(MeV/c 2 ) 


AlnrriH = +0.693 


-41.3 


Am, = +1.8 GeV/c 2 [9] 


11.0 


Aoiem{Q = m z c 2 ) = +0.00035 [10] 


-6.2 


Am z = +2.1 MeV/c 2 [11] 


2.6 



TABLE I: The effect on mw of +1<t increases of the in- 
put parameters dominating the uncertainty on the mw 
prediction. Since the Higgs boson has not been observed, 
we show the effect of doubling the Higgs boson mass from 
100 GeV/c 2 to 200 GeV/c 2 [|. 



Experiment mw (GeV/c 2 ) 



ALEPH [12] 


80.440 


± 


0.051 


OPAL [13] 


80.416 


± 


0.053 


L3 [14] 


80.270 


± 


0.055 


DELPHI [15] 


80.336 


± 


0.067 


CDF Run I [16] 


80.433 


± 


0.079 


D0 Run I [17, .18] 


80.483 


± 


0.084 


LEP Average [19] 


80.376 


± 


0.033 


Tevatron Run I Average 


[20] 80.456 


± 


0.059 


World Average 


80.392 


± 


0.029 



TABLE II: Direct measurements of the W boson mass, 
the preliminary combined LEP average, the combined 
Tevatron Run I average, and the preliminary world av- 
erage. 

The leptonic decays of singly produced Z bosons pro- 
vide important control samples, since both leptons 
from Z boson decay are well measured. The produc- 
tion and decay uncertainties on the measurement of 
mw from pp and e + e~ collider data are almost com- 
pletely independent [2^ ]. 

We present in this Section an overview of W and 
Z boson production at the Tevatron, a description 
of the coordinate definitions and symbol conventions 
used for this measurement, and a broad discussion of 
our mw measurement strategy. 



A. W and Z Boson Production and Decay 

W and Z bosons are produced in ^fs = 1.96 TeV 
pp collisions primarily through s— channel annihila- 
tion of valence u and/or d quarks (Fig. 2]), with a 
smaller O(20%) contribution from sea quarks. The 
quark (antiquark) has a fraction x p (xp) of the pro- 



6 



ton's (antiproton's) total momentum, producing a W 
or Z boson at center of mass energy \fi = Q equal 
to its mass times c 2 . The rate of production can be 
predicted from two components: (1) the momentum 
fraction distributions of the quarks, f a (x, Q 2 ), which 
are determined from fits to world data [23|, [24[ ; and 
(2) a perturbative calculation of the qq' — > W or Z 
boson process (251 ] . 

u 



u 




d(u) 



FIG. 4: Leading-order annihilation of a quark and an- 
tiquark inside the proton and antiproton, respectively, 
producing a W + or Z° boson. The quark (antiquark) 
has energy x p E p (x p Ep), where E p (E p ) represents the 
total proton (antiproton) energy. The production occurs 
at a partonic center-of-mass energy Q. The uu — ► Z° and 
du — > W~ processes are similar. 



W and Z bosons can decay to lepton or quark 
pairs. Decays to quark pairs are not observable 
given the large direct qq' background, and decays to 
t — > !/,-+hadrons arc not as precisely measured as 
boson decays to electrons or muons. For these rea- 
sons we restrict ourselves to the direct electronic and 
muonic decays (W — > ev, W — > [lis, Z — > ee, and 
Z — > fXfi), with the corresponding decays to r — > lep- 
tons considered as backgrounds to these processes 
(Section IVIIip . The branching ratio for each lep- 
tonic decay W -> Iv (Z -> JZ) is (3.3%), and 

the measured cross section times branching ratio is 
(2749 ± 174) pb [(254.9 ± 16.2) pb] [26|. 

B. Conventions 

We use both Cartesian and cylindrical coordinate 
systems, in which +z points in the direction of the 
proton beam (east) and the origin is at the center 
of the detector. In the right-handed Cartesian coor- 
dinate system, +x points north (outward from the 
ring) and +y points upwards; in the cylindrical sys- 
tem, </> is the azimuthal angle and r is the radius from 
the center of the detector in the x — y plane. The 
rapidity y = — |ln[(E — p z c)/(E + p z c)] is additive 




FIG. 5: A IT boson event, with the recoil hadron mo- 
mentum (ut) separated into axes parallel («||) and per- 
pendicular (u±) to the charged lepton. 



under Lorentz boosts along the z axis. For massless 
particles, this quantity is equal to the pseudorapidity 
77 = — ln[tan(0/2)], where 8 is the polar angle with 
respect to the z axis. All angles are quoted in radians 
unless otherwise indicated. 

Because the interacting quarks' longitudinal mo- 
menta p z are not known for each event, we gener- 
ally work with momenta transverse to the beam line. 
The interacting protons and antiprotons have no net 
transverse momentum. Electron energy (muon mo- 
mentum) measured using the calorimeter (tracker) is 
denoted as E (p), and the corresponding transverse 
momenta pr are derived using the measured track 
direction and neglecting particle masses. The event 
calorimetric pr, excluding the lepton(s), is calculated 
assuming massless particles using calorimeter tower 
energies (Section IIII A2[) and the lepton production 
vertex, and provides a measurement of the recoil mo- 
mentum vector ut- The component of recoil pro- 
jected along the lepton direction is denoted un and 
the orthogonal component is u± (Fig. [5]) . The trans- 
verse momentum imbalance in a IT boson event is 
a measure of the neutrino transverse momentum p^ 
and is given by p — ~(p T + Ut), where p^ is the 
measured charged lepton transverse momentum. 

When electromagnetic charge is not indicated, 
both charges are considered. We use units where 
H = c = 1 for the remainder of this paper. 

C. Measurement Strategy 

The measurement of the final state from W — > Iv 
decays involves a measurement of p^ and the total 
recoil Ut- The neutrino escapes detection and the 
unknown initial partonic p z precludes the use of p z 
conservation in the measurement. The boson invari- 
ant mass is thus not reconstructable; rather, the 2- 
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dimensional "transverse mass" m-r is used in the mw 
fit: 



rn,T = y2p^ft T (1 - cos A<j>), (3) 

where A(f> is the angle in the transverse plane between 
the leptons, whose masses are negligible. The fit to 
the rriT distribution provides the statistically most 
precise measurement of mw ■ 

The charged lepton, which can be measured pre- 
cisely, carries most of the observable mass informa- 
tion in the event. We calibrate the muon momen- 
tum using high statistics samples of the meson de- 
cays J ftp — ► and T — > which are fully re- 
constructable and have well known masses. This re- 
sults in a precise track momentum calibration, which 
we transfer to the calorimeter with a fit to the ra- 
tio of calorimeter energy to track momentum (E/p) 
of electrons from W boson decays. The accuracy of 
these calibrations is demonstrated by applying them 
to measurements of the Z boson mass in the muon 
and electron decay channels. We then incorporate 
the known Z boson mass as an additional calibration 
constraint. 

The other directly measurable quantity needed for 
the calculation of tut is the recoil transverse momen- 
tum ut- Since the W and Z bosons are produced at a 
similar Q 2 , they have similar recoil distributions. We 
use the leptons from the Z boson decay to measure 
the pt of the Z boson. We then calibrate our model 
of ut by measuring the balance between the recoil 
and Z boson pr- The Z boson statistics are suffi- 
cient to perform a recoil calibration to 1% accuracy, 
which leads to a systematic uncertainty commensu- 
rate with other uncertainties on mw- 

To accurately model the shape of the tot distri- 
bution, we use a fast Monte Carlo simulation of the 
pp — > W — ► Iv process including the recoil and the 
detector response. The custom fast simulation allows 
flexibility in parametrizing the detector response and 
in separating the effects of the detector model com- 
ponents. We use a binned likelihood to fit the mea- 
sured tot distributions to templates (Section III Dp 
generated from the fast simulation, with mw as the 
free parameter. All mw and lepton energy scale fits 
are performed with this procedure. 

Though less statistically precise, the p l T and 
p T distributions provide additional information on 
the W boson mass and are used as important tests 
of consistency. We separately fit these distributions 
for mw and combine all fits in our final result. 

During the measurement process, all W boson 
mass fits were offset by a single unknown random 
number chosen from a flat distribution in the range 
[-100,100] MeV. The fit result was thus blinded to 



the authors until the analysis was complete [27J • The 
final measured mw and its uncertainty have not 
changed since the random offset was removed from 
the fit results. 

We give a brief overview of the template likelihood 
fitting procedure in Sectkm UlDI Scction HTll describes 
the detector and the fast detector simulation used in 
the analysis. The W boson measurement samples 
are defined in Section [IV] We describe the precision 
measurements of muons and electrons in Sections IVl 
and I VII respectively. These sections include event se- 
lection, calibration, and resolution studies from the 
dilepton and W boson data samples. Measurement of 
the recoil response and resolution is presented in Sec- 
tion IVIII The backgrounds to the W boson sample 
are discussed in Section [Villi Theoretical aspects of 
W and Z boson production and decay, including con- 
straints from the current data sample, are described 
in Section Hxl We present the W boson mass fits and 
cross-checks in Section [X] Finally, in Section IXII we 
show the result of combining our measurement with 
previous measurements, and the corresponding impli- 
cations on the predicted standard model Higgs boson 
mass. 



D. Template Likelihood Fits 

All the fits involving mass measurements and the 
energy scale (Sections El EH and [Xj) are performed 
with a template binned likelihood fitting procedure. 
A given distribution to be fit is generated as a discrete 
function of the fit parameter, using the fast simula- 
tion. These simulated distributions are referred to 
as "templates." For each value of the fit parameter, 
the simulated distribution is compared to the data 
distribution and the logarithm of a binned likelihood 
is calculated. The binned likelihood is the Poisson 
probability for each bin to contain the rij observed 
data events given to,- expected events, multiplied over 
the N bins in the fit range: 

iV _ m rrii 
i—l 

We calculate the logarithm of the likelihood using the 
approximation Inn! w (n + 1/2) ln(n + 1) — n: 

N 

lnX « InTOj — mi — (rt,; + 1/2) ln(n.j + 1) + rij]. 

i=l 

(5) 

The best-fit value of the parameter maximizes the 
likelihood (or equivalently minimizes — In/!), and the 
±l(j values arc those that increase — ln£ by 1/2. The 
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approximation for In nl only affects the shape of the 
likelihood about the minimum and not the position of 
the minimum. The procedure is validated by fitting 
simulated data ( "pseudoexperiments" ) and no bias is 
found. We symmetrize the uncertainty by taking half 
the difference between the +la and —la values. For 
the E/p fits in the W boson sample, we reduce the 
effect of finite template statistics by fitting — In C to 
a parabola, and extracting the best-fit value and the 
uncertainty from this parabola. 



III. DETECTOR AND MODEL 

The CDF II detector [H, 11] is well suited for the 
m\v measurement. Its high-resolution tracker and 
calorimeter measure individual charged lepton mo- 
menta from W and Z boson decays with a resolution 
of pa 2%. It has similar acceptance and resolution for 
central electrons and muons, giving the two channels 
similar weight in a combined mass measurement. 



A. Detector Components 

The CDF II detector (Fig. [6]) is a multi-purpose de- 
tector consisting of: an inner silicon tracker designed 
to measure the production vertex of charged particles 
with high precision; an outer tracking drift cham- 
ber to measure charged particle momenta; a solenoid 
to provide a uniform 1.4 T magnetic field inside 
the trackers; electromagnetic calorimeters to contain 
and measure electron and photon showers; hadronic 
calorimeters for hadron energy measurements; and a 
muon system to detect muons escaping the calorime- 
ters. The detector information is read out on-line 
and saved for later analysis when event topologies 
consistent with a particular physics process (or class 
of processes) are selected. The read-out decision is 
made with a fast three-level trigger system that has 
high efficiency for selecting the W and Z bosons to 
be used in the offline analysis. 



1. Tracking System 

The silicon tracker (Fig. [7]) consists of three sepa- 
rate detectors: Layer 00, SVX II, and ISL. Layer 00 is 
a single layer of 300 ^m thick sensors attached to the 
beam pipe at a radius of 1.3 cm. Five additional lay- 
ers of sensors at radii ranging from 2.5 cm to 10.6 cm 
comprise SVX II. Surrounding these sensors are port 
cards, which transport deposited charge information 



from the silicon wafers to the readout system. The in- 
termediate silicon layers (ISL) are located at radii of 
20.2 cm and 29.1 cm. The SVX II is segmented longi- 
tudinally into three barrels in the region \z\ < 45 cm. 
This covers the pp interaction region, which is well ap- 
proximated by a Gaussian distribution with a z ~ 30 
cm. We do not use the silicon measurements in this 
analysis, though we model the tracker's effects on lep- 
tons and photons (Section IIII Bp . 

An open-cell drift chamber, the central outer 
tracker (COT) [29| . surrounds the silicon tracker 
and covers the region \z\ < 155 cm (\rj\ < 1) and 
40 cm < r < 137 cm. The COT consists of eight con- 
centric 'super layers,' separated azimuthally into cells. 
Each cell contains 12 sense wires to measure the ion- 
ization produced by a charged particle in the ambient 
argon-ethane gas mixture. The superlayers alternate 
between a purely axial configuration, with sense wires 
parallel to the beam line, and a small-angle stereo 
configuration, with sense wires at a 2° angle relative 
to the z axis. 

The sense wires are strung from end to end in z 
and held under tension at each aluminum endplate 
(Fig. [5]). The wires are azimuthally sandwiched by 
field sheets, which provide a 1.9 kV/cm electric field. 
All cells are rotated at a 35° angle relative to a radial 
line, such that the ionized electrons travel approxi- 
mately azimuthally to the wire under the combined 
influence of the local electric field and the global mag- 
netic field from the solenoid. 

Within a given cell the sense wires are slightly off- 
center relative to the field sheets. In addition, the 
sense wires and field sheets sag under the influence 
of gravity, with the field sheets sagging more due to 
their larger masses. These effects cause a small elec- 
trostatic deflection of the sense wires toward a par- 
ticular field sheet. To prevent the relative deflection 
of sense wires within a cell, a support rod connects 
the sense wires at the center of the detector. The 
support rod results in a small (« 2 mm) region at 
z = cm where charged particles are not measured. 

Between the solenoid and the COT is a time-of- 
flight system (TOF) consisting of scintillator bars 
that precisely measure the time of incidence of 
charged particles. From this measurement and the 
tracker information, a particle's velocity and mass 
can be inferred. The TOF is not utilized in this anal- 
ysis. 



2. Calorimeter System 

The CDF calorimeter is segmented radially into 
electromagnetic and hadronic sections. The central 
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calorimeter covers \rj\ < 1.1 and is split at the cen- 
ter into two separate barrels covering +77 and —rj. 
Each barrel consists of 24 azimuthal "wedges" of 
size 0.26 radians (15°) with ten projective towers of 
size A77 w 0.11. To allow a pathway for the solenoid 
cryogenic tubes, a two-tower region is removed, cor- 
responding to 0.77 < 77 < 1.0, 75° < <j) < 90°, and 
z > 193 cm. The forward calorimeter covers 
1.1 < 1 77 I < 3.6, filling the forward gaps with a plug 
shape (Fig.jnj). 

The central electromagnetic calorimeter (CEM) 
[3C| has a thickness of 18 radiation lengths, consist- 
ing of 31 radial scintillator layers interleaved with 30 
layers of lead-aluminum plates. At a radius of 184 
cm electromagnetic showers have traversed about six 
radiation lengths (including the solenoidal coil) and 
have their maximum energy deposition. At this ra- 
dius finely segmented strip and wire chambers (CES) 
measure the energy deposition with a position reso- 
lution of Rj 2 mm. 

The local shower position in the azimuthal direc- 
tion in the tower is denoted as CES x, which ranges 



from -24.1 cm to 24.1 cm. The wire chambers extend 
only to \x\ < 22.5 cm, and for \x\ > 23.1 cm no energy 
measurements are made. In this region wavelength 
shifters read out the light from the scintillator, and 
steel and foam separate the towers. Light follows a 
waveguide to a phototube positioned at the back end 
of the hadronic calorimeter. 

Parallel to the beam line, the position at shower 
maximum is denoted CES z. The strip chambers ex- 
tend from 6-239 cm in \z\, and there is no scintillator 
for \z\ < 4.2 cm, where the two calorimeter barrels 
meet . 

The central hadronic calorimeter [3l[ is separated 
into a central region (CHA, |?7| < 0.6) with 32 longitu- 
dinal layers of scintillator sandwiched with steel and 
a forward "wall" calorimeter (WHA, 0.6 < \t)\ < 1.1) 
with 15 such layers. These calorimeters have thick- 
nesses of « 4.5 interaction lengths. 

The plug calorimeter [32| has a comparable de- 
sign to the central calorimeter with scintillator- 
lead electromagnetic calorimeters and scintillator- 
steel hadronic calorimeter compartments. The <j> seg- 
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FIG. 7: End view of the silicon detector. The inner- 
most layer (Layer 00) is attached to the beam pipe, and 
is surrounded by five concentric layers of silicon wafers 
(SVX II). The outermost layers are the intermediate sil- 
icon layers (ISL) , which sit just inside the outer tracking 
chamber. 



mentation is 0.13 radians up to \rj\ = 2.1, and then 
broadens to 0.26 radians. The two furthest forward 
plug towers cover the regions 2.6—3.0 and 3.0—3.6, 
while the remaining towers have a size Arj = 0.1. 



3. Muon Detectors 

The muon systems relevant for the W mass mea- 
surement cover the region < 1. The central 
muon detector (CMU) and the central muon upgrade 
(CMP) cover \rj\ < 0.6, while the central muon exten- 
sion (CMX) covers 0.6 < |ry| < 1. 

The CMU detector 33] is located at the outer edge 
of the CHA, 347 cm from the z axis. The CMU is 
segmented into 15° azimuthal wedges containing four 
layers of proportional drift chambers that cover 12.6°. 
The maximum drift time within a chamber is 800 ns, 
about twice as long as the 396 ns spacing between 
pp crossings. CMU information must therefore be 
combined with reconstructed COT particle tracks to 
determine the appropriate pp crossing. 

Because the total thickness of the central calorime- 
ter is about five interaction lengths, approximately 
0.5% of high-momentum pions reach the CMU. To 
reduce this background, the CMP detector is located 
behind an additional 60 cm of steel. The CMP has 
a similar construction to the CMU, with the excep- 
tion that wider drift chambers are used to cover the 
same solid angle, resulting in a maximum drift time 



of 1.8 /us rather than 800 ns. 

The CMX detector ]34] consists of eight drift 
chamber layers beyond both the calorimeter and the 
steel detector support structure (6 — 10 interaction 
lengths). The CMX <f> regions used in this analysis 
are -45° < cj) < 75° and 105° < < 225°. New 
detectors for Run II cover much of the remaining <f> 
region, but were not fully commissioned for the data- 
taking period of this analysis. Scintillator detectors 
(CSX) at the inner and outer surfaces of the CMX 
provide timing information to the trigger to separate 
collision particles from other sources such as beam 
halo or cosmic rays. 



4- Trigger System 

The trigger consists of three stages with progres- 
sively greater sophistication of event reconstruction. 
The first stage is hardware-based, the second a mix 
of hardware and software, and the third a farm of 
processors performing full event reconstruction. 

The first trigger stage, level 1, includes tracker, 
calorimeter, and muon reconstruction. The charged 
particle track reconstruction is performed with the 
extremely fast tracker (XFT) [35| based on the four 
axial COT superlayers. A track segment is recon- 
structed in a given superlayer if at least 11 of the 12 
sense wires [36[ in a wide road have charge deposi- 
tion above a given threshold ( "hits" ) . The list of seg- 
ments from the full tracker is compared to predefined 
groups of segments expected from charged particles 
above a given momentum threshold. When matches 
are found, track candidates are created and passed 
to the track extrapolator (XTRP) [H . The XTRP 
determines the expected <f> positions of the tracks in 
the calorimeter and muon detectors, for the purpose 
of forming electron and muon candidates. 

The calorimeter reconstruction at level 1 defines 
separate electromagnetic and hadronic "trigger tow- 
ers" as tower pairs adjacent in r\. The tower px is 
calculated assuming a collision vertex z = and an 
electron candidate is formed if the ratio of hadronic to 
electromagnetic energy (Had/EM) in a trigger tower 
is less than 1/8. The high-momentum electron trigger 
used in this analysis requires a level 1 trigger tower 
with electromagnetic px > 8 GeV matched to a track 
with pt > 8 GeV, and drops the Had/EM require- 
ment for electromagnetic px > 14 GeV. 

Level 1 muon reconstruction includes a px esti- 
mate within the CMU and CMX chambers from 
the relative timing of the hits in different layers. 
The CMU track segments are combined with recon- 
structed CMP track segments to create "CMUP" 
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FIG. 8: End view of a section of a central outer tracker (COT) endplate. The COT consists of eight concentric 
'superlayers,' separated azimuthally into cells, each containing 12 sense wires and sandwiched by field sheets. The 
endplates contain precision-machined slots where each cell's sense wires and field sheets are held under tension. The 
radius at the center of each superlayer is shown in cm. 



muon candidates. For the majority of the data CMX 
candidates also require local CSX hits consistent with 
particles originating from the collision. For our W 
and Z boson samples we use a muon trigger that re- 
quires CMU or CMX px > 6 GeV matched to an 
XFT track with p T > 4 GeV (CMUP) or p T > 8 
GeV (CMX). 

The level 2 calorimeter reconstruction uses a more 
sophisticated clustering algorithm for electromag- 
netic objects. This improves energy measurement 
resolution and allows a higher threshold (pt > 16 
GeV) to be applied. To reduce rates, the XFT 
track requirement for CMUP candidates was raised 
to px > 8 GeV for most of the data-taking period. 

At level 3, approximately 300 dual processor com- 
puters allow full track pattern recognition, muon re- 
construction, and calorimeter clustering. Variables 
used to select electrons at level 3 are the lateral 
shower profile, L s f lr (Section IIV Bp . and the distance 
between CES z and the z-position of the track extrap- 
olated to the CES (Az). The L s h r variable quantifies 
the difference between the measured energies of tow- 
ers neighboring the electron in rj and the expected 
energies determined from electron test beam data. 
The trigger requirements of L s h r < 0.4 and |Az| < 8 
cm are ~100% efficient for electrons from W and Z 
boson decays. The high-momentum electron trigger 



also requires electromagnetic px > 18 GeV and track 
Pt > 9 GeV. For efficiency studies we use a separate 
trigger that requires electromagnetic pt > 25 GeV 
and p^ 3 > 25 GeV, but has no quality requirements 
at level 3 and no trigger track requirements. At level 
3, is defined as the negative of the vector sum 
of the transverse momenta in all calorimeter towers. 
The high-momentum muon trigger requires a COT 
track with p T > 18 GeV matched to a CMUP or 
CMX track segment. 



5. Luminosity Detector 

The small-angle Cherenkov luminosity counters 
(CLC) [H[ are used to measure the instantaneous 
and integrated luminosity of our data samples. The 
CLC consists of two modules installed around the 
beampipe at each end of the detector, providing cov- 
erage in the regions 3.6 < \rj\ < 4.6. Each module 
contains 48 conical gas Cherenkov counters pointing 
to the collision region. Signals in both CLC modules 
coinciding in time with the bunch crossing are used to 
measure the instantaneous luminosity and to trigger 
collision events. Events collected with this trigger, 
known as "minimum bias" events, are used to study 
the detector response to generic inelastic pp collisions 
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(Section En). 



Ionization Energy Loss 



B. Detector Model 



We use a parametrized model of the detector re- 
sponse to electrons, muons, and the hadronic recoil. 
The model is incorporated into a custom fast simu- 
lation that includes lepton and recoil reconstruction, 
event selection, and fit template generation. The sim- 
ulation provides both flexibility in determining the ef- 
fects of various inputs, and computing speed to allow 
frequent high-statistics studies. A sample of O(10 7 ) 
events can be generated using a single-processor ma- 
chine in one day. This is several orders of magnitude 
more than the O(10 3 ) events that can be produced 
with the standard GEANT-based CDF simulation [401 

B 

We describe in this section the simulation of elec- 
trons and muons. Fits to the data that determine the 
values of some of the model parameters are described 
in Sees. [V] and I VII The detector model of hadronic 
recoil response and resolution is discussed in Sec. lVIH 

The model components common to muons and 
electrons are: ionization energy loss and multiple 
scattering in the beam pipe and tracker volume; 
parametrized track hit resolutions and efficiencies; 
and track reconstruction. We describe these com- 
ponents in the muon simulation overview, and then 
discuss the electron- and photon-specific simulation. 



Muon Simulation 



The differential ionization energy loss of muons and 
electrons in the tracking system is simulated accord- 
ing to the Bethe-Bloch equation [Til ]: 



dE 
dx 



KZ 
Aft 



1 2?n e /3 T max 

2 n (1-/3 2 )/ 2 



a 2 



(6) 



where K = AirN AfeTn e , Na is Avogadro's number, 
r e is the classical electron radius, Z(A) is the atomic 
(mass) number, (3 is the particle velocity, I is the 
mean excitation energy, T max is the maximum kine- 
matic energy that can be transferred to a free electron 
in a single collision, and 8 is the material-dependent 
density effect as a function of (3 [ll[ . When calculat- 
ing the effect of S, we take the material to be silicon 
throughout. 

To calculate muon energy loss in the material up- 
stream of the COT (r < 40 cm), we use a three- 
dimensional lookup table of the material properties 
of the beam pipe, the silicon detector, and the wall 
of the alumnium can at the inner radius of the COT. 
The lookup table determines the appropriate Z/A 
and / values, along with the radiation length X ( Ap- 
pendix [A"| . for each of 32 radial layers. Except for the 
inner and outer layers, the map is finely segmented 
longitudinally and in azimuth to cap ture the material 
variation in the silicon detector (42J • Inside the COT 
fiducial volume we calculate the energy loss between 
each of the 96 radial sense wires. 

The energy loss model is tuned using the data. We 
apply a global correction factor of 0.94 to the calcu- 
lated energy loss in the material upstream of the COT 
in order to obtain a J/ip — > fj,fj, mass measurement 
that is independent of the mean inverse momentum 
of the decay muons (Section IV B 3[) . 



Muon and electron tracks are reconstructed using 
only COT hit and beam position information (Sec- 
tion [IV]). Thus, the simulation of the silicon detector 
consists entirely of energy loss and multiple scatter- 
ing. In the COT, hit resolutions and efficiencies are 
additionally simulated, and track reconstruction is 
performed. The total measured muon EM calorime- 
ter energy is simulated by combining the minimum- 
ionizing energy deposition with energy from final- 
state photon radiation (Section IIXD[) and the recoil 
and underlying event [39|]. Finally, the detector fidu- 
ciality of muons is calculated using a map of the muon 
detector geometry as a function of r\ and <f>. The map 
is extracted from a full GEANT-based simulation of 
the CDF II detector SSI. 



Multiple Coulomb Scattering 

Multiple Coulomb scattering in the beampipe, sil- 
icon detector, and COT affects the resolution of the 
reconstructed track parameters for low-momentum 
tracks. We model the scattering u sing a Gaussian 
distribution for 98% of the scatters [43[ with an an- 
gular resolution oq defined by 



13.6 MeV 

Wp 



(7) 



where x is the thickness of the layer and Ao is the 
layer's radiation length (Section IIIIB2p . Simulation 
of multiple scattering is implemented for each radial 
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layer of the three-dimensional lookup table and be- 
tween each COT layer. 

Based on the results of low-energy muon scatter- 
ing data |44|], we model the non-Gaussian wide-angle 
scatters by increasing ag by a factor of 3.8 for 2% of 
the scatters. 



COT Simulation and Reconstruction 

The charged track measurement is modeled with 
a full hit-level simulation of the charge deposition in 
the COT and a helical track fit. The parameter reso- 
lution of reconstructed tracks is affected by the indi- 
vidual hit resolution, and by the distribution of the 
number of hits (Nhu) used in the fit [45]. 

We tune the COT hit resolution using the width of 
the T — * fifi mass distribution reconstructed with 
non-beam-constrained tracks. The tuned value of 
[150±3(stat)] (xm is consistent with the 149 /im RMS 
of the observed hit residual distribution for the muon 
tracks in Z — > fj.fi data. We use a 150 ^m hit resolu- 
tion for the simulation of the T, W, and Z bosons. 

We use a dual-resolution model to describe the nar- 
rower mass peak in the high-statistics J/ip — > fifi 
sample, where the muons generally have lower mo- 
menta than the other samples. The J/ip mass peak 
width is particularly sensitive to multiple scattering 
and relative energy loss, and our hit-resolution model 
compensates for any mismodcling that affects the 
peak width. We find that a single-hit resolution of 
155 fim applied to 70% of the tracks and 175 /xm ap- 
plied to the remaining 30% adequately describes the 
width and lineshape of the J/ip —* fx\i mass peak. 

To describe the N} lit distribution, we use a dual- 
hit-efficiency model, the larger one applied to the ma- 
jority of the tracks. The lower efficiency accounts for 
events with high COT occupancy, where fewer hits 
are attached to reconstructed tracks. The two pa- 
rameters are tuned to match the mean and RMS of 
the data Nhu distributions. We independently tune 
these parameters for the J/ip sample, the T sample, 
and the W and Z boson samples. 

COT hit positions from a charged track are used 
to reconstruct a helix with a ^-minimization proce- 
dure. The axial helix parameters [46[ are the impact 
parameter with respect to the nominal beam position, 
do, the azimuthal angle at the closest approach to the 
beam, </>o, and the curvature of the track, c, defined 
to be (2R)~ 1 , where R is the radius of curvature. 
The stereo helix parameters are the longitudinal po- 
sition at the closest approach to the beam, zq, and 
the cotangent of the polar angle, cot 0. 

When optimizing resolution of lepton tracks from 



prompt resonance decays, we constrain the helix to 
originate from the location of the beam. The trans- 
verse size of the beam is « 30 fxrn at z — cm 
and increases to 50 — 60 /j,m at \z\ = 40 cm [47]]. 
For simplicity we assume an average beam size of 
[39 ± 3(stat)] /im, which is determined from a fit 
to the width of the Z — > fifi mass peak. The 
beam constraint improves the intrinsic fractional mo- 
mentum resolution by about a factor of three, to 
Spt/pt ~ 0.0005 p T /GeV. 

We perform a track fit on our simulated hits in the 
same manner as the data. The hits are first fit to a 
helix without a beam constraint; hits with large resid- 
uals (> 600 fim) are dropped from the track (in order 
to remove spurious hits added in data pattern recog- 
nition); and the track is fit again with an optional 
beam constraint. This option is applied to prompt 
lepton tracks from W and Z boson decays, but not to 
tracks from J/ip decays, approximately 20% of which 
are not prompt. The prompt muons from T decays 
are fit twice, both with and without the beam con- 
straint, as a consistency check. 



Calorimeter Response 

Muons deposit ionization energy in the calorimeter. 
We simulate a muon's EM energy deposition using 
a distribution taken from cosmic ray muons passing 
through the center of the detector, in events with 
no other track activity. An additional contribution 
comes from energ y flo w into the calorimeter from the 
underlying event [391 ] . We model this energy using a 
distribution taken from W — > fiv data events, using 
towers separated in azimuth from the muon. 

Muons with a CES z position within 1.58 cm of 
a tower boundary typically deposit energy in two 
calorimeter towers. We use this criterion in the sim- 
ulation to apply the underlying event and final-state 
photon radiation ([Section llXDp contributions for one 
or two towers. The simulated underlying event en- 
ergy includes its dependence on Hi and u± (Fig. O, 
and on the tower 77 position of the muon when it 
crosses the CES (Section EHH). 



Detector Fiduciality 

The CMUP and CMX muon systems do not have 
complete azimuthal or polar angle coverage. We cre- 
ate an rj — <p map of each muon detector's coverage 
using muons simulated p~fl ] with a detector geometry 
based on geant [4(|. We use the map in the fast 
simulation to determine the fiduciality of a muon at 
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a given r\ — qb position. 

We incorporate the relative efficiency of the CMUP 
to CMX triggers in the fast simulation by matching 
the ratio of CMUP to CMX events in the W -> fiv 
data (Section ItTA"]) . 



2. Electron and Photon Simulation 

The dominant calibration of the calorimeter energy 
measurement E of electrons uses their track momenta 
p and a fit to the peak of the E/p distribution. An 
additional calibration results from a mass fit to the 
Z boson resonance and reduces the calibration uncer- 
tainty by 20% relative to the E/p calibration alone. 

The E/p method relies on an accurate modeling 
of radiative effects that reduce the track momentum 
measured in the COT. A given electron loses « 20% 
of its energy through bremsstrahlung radiation in the 
silicon detector, and this process has the most sig- 
nificant impact on the E/p calibration. The total 
amount of silicon detector material is tuned with data 
using highly radiative electrons (Section IVI A[) . We 
additionally model processes that affect the shape 
of the E/p distribution: photon conversion in the 
tracker; energy loss in the solenoid and the time-of- 
flight system; electromagnetic calorimeter response 
and resolution; and energy loss into the hadronic 
calorimeter. The models of ionization energy loss and 
multiple scattering in the tracker, as well as the COT 
track simulation and reconstruction, are the same as 
for muons (Section IIII B ip . 



Bremsstrahlung 



The differential cross section for an electron of en- 
ergy E e to radiate a photon of energy E~, is given by 
the screened Bethe-Heitler equation [48j over most of 
the y = E-y/E e spectrum. In terms of the material's 
radiation length Xo, the differential cross section for 
bremsstrahlung radiation is: 



da A 
dy NaXqp 



y 



(8) 



where C is a small material-dependent correction 
(Appendix [A"|) . Figure [5] shows the integrated thick- 
ness of material upstream of the COT, in terms of ra- 
diation lengths, traversed by the reconstructed elec- 
tron tracks in W — > eu data. The number of photons 
emitted per layer is given by: 



n t = JT 



l + c) (y - lny -l) + \{l~y ) 2 
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FIG. 9: The distribution of material upstream of the COT 
traversed by reconstructed electron trajectories in W — > 
eu data events, in units of radiation lengths. The peaks 
at « 0.08, w 0.13, and « 0.24 correspond respectively 
to trajectories outside the silicon detector {\z\ > 45 cm), 
within the silicon detector, and crossing silicon barrels 
(\z\ S3 15,45 cm). The mean of the distribution is 19%. 



where x is the thickness of the layer and yo is a lower 
threshold introduced to avoid infrared divergences. 
We use y = 10" 4 [I9J and determine C = 0.0253 
using the silicon atomic number Z = 14. 

For each layer of the silicon or COT material, we 
use a Poisson distribution with mean iV 7 to determine 
the number of photons radiated in that layer. For 
each radiated photon, we calculate y from the spec- 
trum in Eq. ([8]). To correct for inaccuracies of the 
screened Bethe-Heitler equation at the ends of the y 
spectrum, we apply a suppression factor if y < 0.005 
or y > 0.8. 

For radiation of high-momentum photons (y > 
0.8), the approximation of complete screening of the 
nuclear electromagnetic field by the atomic electrons 
breaks down. In this region, the full Bethe-Heitler 
equation for incomplete screening [48| must be used. 
We implement this correction by removing generated 
photons in the high-?/ region such that we match the 
reduced cross section from incomplete screening. 

Two effects reduce the cross section for low- 
momentum photon radiation 50] : multiple scattering 
and Compton scattering. Multiple Coulomb scatter- 
ing suppresses long-distance interactions, and the re- 
sulting LPM suppression (5lj in low-momentum ra- 
diation can be expressed in terms of the Bethe-Heitler 
cross section 521: 



(9) 
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where Elpm depends on the material. We use 
Elpm = 72 TeV, appropriate for silicon, and apply 
the suppression when Slpm < 1- 

Radiated photons scatter off the atomic electrons, 
and destructive interference of low-momentum pho- 
tons suppresses this radiation (53|. The suppression 
factor is: 



S, 



C ompton 



v 



Ei/Er 



(11) 



where E p — jutp is 2.4 MeV for a 40 GeV electron in 
silicon, using the silicon plasma frequency u> p , and 7 
is the Lorentz factor. 

In any given simulated event, the product of Slpm 
and Scompton provides the probability that a pho- 
ton generated from the screened Bethe-Heitler equa- 
tion with y < 0.005 survives the low-momentum 
suppression. For a 40 GeV electron radiating a 20 
MeV (8 MeV) photon, the suppression factors are 
Slpm = 0.95 (0.60) and S Compton = 0.99 (0.92). 
Our simulated y spectrum from W boson decay elec- 
trons reproduces the spectrum obtained by a GEANT 
14011 simulation. 



Photon Conversion 

Photons can convert to an electron-positron pair by 
interacting with the tracker material. The differential 
cross section for a photon of energy E 1 > 1 GeV to 
convert into an electron with energy E e is given by 
the screened Bethe-Heitler equation |48j : 



da _ A 
dy N A X a p 



[l-(4/3 + C)y(l-y)] 



(12) 



where y = E e /E 1 . Integrating over y and multiplying 
by pxNa/A gives the total cross section, from which 
we obtain the following conversion probability at high 
photon energy: 



P, 



s - (Ej — ► 00) = 1 — e 



-(7/9-C/6)x/X 



(13) 



We parametrize the cross section as a function of pho- 
ton energy using the tables for photon cross sections 
in silicon given in [54] . We apply the ratio shown in 
Fig. [10] to the high-energy cross section when calcu- 
lating the conversion probability. 

For each radiated photon upstream of the COT, 
we integrate the material between the radiation point 
and the COT inner can. If the photon converts, we 
take the conversion point to be halfway between the 
radiation point and the inner can. If the photon does 
not convert before the COT, we integrate the ma- 
terial in the COT and take a converting photon to 
convert halfway through the COT. 
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FIG. 10: The ratio of the photon conversion cross section 
at a given photon energy to the cross section at i? 7 = 
100 GeV [5J]. We use this function to scale down the 
cross section obtained from the Bethe-Heitler equation 
[Eq. (fTJ] HI. 



We use the conversion electron momentum spec- 
trum from Eq. (|12[) . ignoring the small effect of the C 
term on the shape. If a radiated photon has high mo- 
mentum, a conversion electron's measured momen- 
tum can be larger than that of the electron from the 
W boson decay. To mimic the offline reconstruction, 
we assign the track from the highest momentum elec- 
tron to the electron cluster. 



Compton Scattering 

The cross section for a low-momentum photon to 
scatter off an electron is similar to that of conversion 
into an e + e~ pair. The differential cross section with 
respect to the photon fractional energy loss y can be 
approximated as (Appendix [A")) : 



dy 



lh 



(14) 



Using a lower bound of y — 0.001, this spectrum 
approximates the Compton energy loss distribution 
for photons radiated from electrons from W boson 
decays. 

We calculate the total cross section in terms of the 
pair production cross section using the tables for pho- 
ton interactions in silicon in [54j . The ratio of cross 
sections as a function of energy is parametrized as 
(Fig. ED : 
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FIG. 11: The ratio of the Compton scattering cross sec- 
tion at a given photon energy to the pair-production cross 
section at E-y = 100 GeV [5J]. This ratio is applicable for 
photons traversing silicon. 



where F{E 1 ) = 2.35e- 116B -< + 2A2e~ lbAE -< - 5.21 - 
0.151-B 7 , with E 1 in GeV, and 7 and 7' are the initial- 
and final-state photons, respectively. We thus use the 
following Compton scattering probability per layer: 

P 7 ^y = R C om(7/9 - C/6)x/X . (16) 



Energy Loss in Solenoid 

After exiting the tracker electrons and photons 
travel through the time-of-flight (TOF) system and 
the solenoid. These systems have thicknesses of 
« 10% and « 85% of a radiation length, respec- 
tively. With this much material it becomes pro- 
hibitive to model individual radiative processes, and 
we instead use a parametrized ener gy-l oss model de- 
termined from a GEANT simulation [4CJ]. The energy 
loss is defined as the difference in energy of a sin- 
gle particle entering the TOF and the total energy of 
particles exiting the solenoid. 

Figure [T^l shows the mean energy loss as a function 
of log 10 (j?T/GeV) of the incoming particle for both 
photons and electrons. Electrons lose more energy 
than photons due to their ionization of the material. 
Since electrons with pr < 400 MeV curve back to the 
center of the detector before exiting the solenoid, we 
do not parametrize energy loss in this energy region. 

The energy loss distribution at a given particle px 
is reasonably described by an exponential. We use 
this distribution, with a mean determined by Fig. [T^l 




log 10 [/» r (GeV)] 



FIG. 12: The mean pr loss as a function of 
log 10 (pT/GeV) for electrons with pr > 400 MeV and pho- 
tons traversing the time-of-flight system and solenoid. 



to model the energy loss of a given particle passing 
through the TOF and solenoid. 



Calorimeter Response and Fiduciality 

The calorimeter simulation models the response of 
the electromagnetic calorimeter as a function of each 
particle's energy and position, and the fraction of 
shower energy leaking into the hadronic calorimeter. 

The electromagnetic calorimeter response, or the 
average measured energy divided by the true particle 
energy entering the calorimeter, can depend on each 
particle's energy. Possible sources of this dependence 
are variations in light yield as a function of calorime- 
ter depth, attenuation in the light guide from the 
scintillator to the phototube, or leakage of showering 
particles into the hadronic calorimeter. The mean 
fractional energy leakage into the hadronic calorime- 
ter for particles exiting the tracker, determined us- 
ing the GEANT calorimeter simulation, is shown as a 
function of log 10 (pr / GeV) in Fig. Q2] 

For a low-pT particle exiting the tracker, the dis- 
tribution of energy loss into the hadronic calorimeter 
is adequately described by an exponential. For high- 
Pt particles (> 10 GeV), the distribution has a peak 
at non-zero values of energy loss. In this energy re- 
gion we model the hadronic energy loss fluctuations 
with the distributions shown in Fig. [T3J Because a 
non-negligible fraction of electrons lose a significant 
amount of energy (5— 10%) in the hadronic calorime- 
ter, it is important to model the energy loss spectrum 
in addition to the mean hadronic energy loss. 
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FIG. 13: The photon and electron pr leakage into the 
hadronic calorimeter. Top: The mean pr leakage as a 
function of log 10 (pT /GeV). Bottom: The distributions of 
Pt leakage for high-pr (> 10 GeV) photons and electrons. 



To correct for any unaccounted dependence of the 
response on incoming particle energy, we use an em- 
pirical model of response that increases linearly with 
particle pr- 



Rem(pt) = S B [l+£(pr/GeV- 39)]. 



(17) 



We determine the slope parameter £ = [6± 7(stat)] x 
10~ 5 using fits to the electron E/p distribution as 
a function of pt in W — > ev and Z — > ee events 
(Section IVI[) . The inclusive E/p distribution from 
W — > ev events is used to calibrate the absolute 
response Se- Since electrons in this sample have a 
mean px of 39 GeV, the fitted values for Se and 
£ are uncorrelated. The parameter £ describes the 
"non-linearity" of the calorimeter response. 

Light attenuation in the scintillator results in non- 
uniform response as a function of distance from the 



wavelength-shifting light guides. The attenuation 
function was measured using test beam data at con- 
struction, and aging effects are measured insitu us- 
ing electrons from W boson decays. The function is 
parametrized as a quadratic function of the CES x 
position within a tower and corresponds to a reduc- 
tion in response of w 10% at the edge of the tower. 
We simulate the light attenuation by reducing the 
energy deposited by each particle according to this 
function, evaluated at the particle's CES x position. 

To improve measurement resolution in data, we 
correct for attenuation effects by applying the inverse 
of the quadratic attenuation function to the measured 
EM energy. We match this procedure in the simula- 
tion. 

The EM calorimeter response drops rapidly as a 
particle crosses the edge of the scintillator and into 
the dead region between towers [30]. We take the 
calorimeter to have zero response for any particle 
with |CES x\ > 23.1 cm or |CES z\ < 4.2 cm. For the 
mw measurement we only use high-energy electrons 
far from the dead regions (Section IIVB|) . 

We apply the following smearing to the calorimeter 
cluster energy: 



a E /E = V0.135 2 /p T + k 2 



(18) 



where the constant term n is determined to be 
[0.89±0.06(stat)±0.13(sys)]% from a fit to the width 
of the electron E/p peak in W boson decays [Ha ]. We 
find further energy smearing is necessary to model 
the multi-particle energy clusters populating the high 
E/p region. When a simulated W or Z decay elec- 
tron radiates in the tracker, we apply an additional 



fractional resolution of 



.3±2.2(stat)]% to each 



radiated particle. This smearing contributes ss 1.3% 
to the effective constant term, and is determined from 
a fit to the width of the Z boson mass peak recon- 
structed from radiative electrons (E/p > 1.06). 

The final contribution to the electron cluster en- 
ergy comes from the underlying event [3^ ] and ad- 
ditional pp interactions. As with muons, we mea- 
sure this energy distribution in W boson data as a 
function of U\\, u±, and the electron tower r\ (Sec- 
tion IVII Bp . These measurements are incorporated 
in the simulation. 



IV. W BOSON SELECTION 

The W boson samples are collected with triggers 
requiring at least one central (|?7| < 1) lepton can- 
didate in the event. A narrow kinematic region is 
defined for W boson selection: 30 GeV < lepton 
p T < 55 GeV; 30 GeV < p T < 55 GeV; 60 GeV 
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< m T {l,p T ) < 100 GeV; and u T < 15 GeV. This 
selection results in low background while retaining 
events with precise my/ information. Additional 
background rejection is achieved through event se- 
lection targeting the removal of Z boson decays to 
leptons. To minimize bias, lepton selection criteria 
are required to have high efficiency or to be explic- 
itly modeled by our fast simulation. 



A. W —* \w Selection 

Muons are identified based on their reconstructed 
COT track quality and production vertex, minimum 
ionizing energy deposited in the calorimeter, and the 
consistency of the track segments reconstructed in 
the muon chambers with the COT tracks. 

All charged lepton candidates from W and Z bo- 
son decay are required to have fully-fiducial cen- 
tral (\zq I < 60 cm) COT tracks with at least 5 hits 
on each of > 3 axial superlayers and > 3 small- 
angle stereo superlayers. For muon candidates we 
remove background from decays of long-lived hadrons 
to muons ("decays in flight") by requiring the track 
impact parameter to be small (|do| < 1 mm) and the 
track fit quality to be good (x 2 /dof < 3). After this 
initial selection, the COT track parameters are up- 
dated with an additional constraint to the transverse 
position of the beam, which has a size of « 30 /zm in 
the luminous region. The beam constraint results in 
a factor of ~ 3 improvement in momentum resolution 
for muons from W boson decays. 

Each muon candidate's COT track is extrapolated 
to the calorimeter and its energy deposition in the 
electromagnetic and hadronic calorimeters is sepa- 
rately measured. Muons near a tower edge in the 
z direction cross two calorimeter towers, and those 
tower energies are combined to determine the muon's 
total energy deposition. We require the muon's elec- 
tromagnetic energy deposition Eem to be less than 2 
GeV and its hadronic energy deposition Eh ad to be 
less than 6 GeV [H- 

All W muon candidates must have a track segment 
in either the CMU and CMP detectors, or the CMX 
detector. COT tracks extrapolated to these detectors 
must have r — 4> positions that match to within 3,5, 
or 6 cm of the CMU, CMP, or CMX track segment 
positions, respectively. 

The Z/"f* — > fifi process presents a significant 
background to the W — > \xv sample. We reduce 
this background by removing events with a second 
opposite-charge muon candidate passing the above 
selection, or passing the following looser set of cri- 
teria: an opposite-charge track with > 10 GeV, 



\do\ < 1 mm, > 2 axial superlayers with > 5 hits, 
and > 2 (1) small-angle stereo superlayers with > 5 
hits for tracks fully (partially) fiducial to the COT; 
Eem < 2 GeV and Eh ad < 6 GeV; and calorime- 
ter isolation < 0.1. Calorimeter isolation is defined 
as the calorimeter pt in an rj — <f> cone of radius 0.4 
surrounding the muon calorimeter towers, divided by 
the muon track px- For events with one identified 
W decay muon and a second muon candidate pass- 
ing the looser criteria, the identified W decay muon 
must also have isolation < 0.1 for the event to be re- 
jected from the W boson sample. The full W boson 
sample, after kinematic selection and Z boson rejec- 
tion, contains 51,128 events in (190.8 ± 11.1) pb _1 of 
data. 

The identification efficiency of muons has a small 
dependence on the recoil in W — > fiv and Z — > fifi 
events, due primarily to the track x 2 an d do require- 
ments. We measure this dependence using Z — > 
events, selected with one muon passing the W muon 
candidate criteria and a second "probe" muon iden- 
tified as a track with pr > 30 GeV. The two muons 
must have opposite charge and reconstruct to an in- 
variant mass in the 81 — 101 GeV range. The frac- 
tion of probe muons passing the additional W muon 
candidate selection criteria is shown in Fig. [14] as a 
function of net recoil energy along the muon direction 
(u||). The observed dependence is parametrized as: 

e = a[l + &(u|| + |u|||)], (19) 

where a is a normalization factor that does not affect 
the my/ measurement and b = [— 1.32± 0.40(stat)] x 
10 -3 . We vary b by ±3u in simulated data and fit 
for my/- Assuming a linear variation of my/ with b, 
we derive uncertainties of Smyy = 1, 6, and 13 MeV 
for the rriT, Pt, and p T fits, respectively. 

B. W —* ev Selection 

Electron identification uses information from the 
COT track quality and production vertex, the match- 
ing of the track to calorimeter energy and position, 
and the longitudinal and lateral calorimeter energy 
profiles. 

An electron candidate's COT track has the same 
fiduciality and hit usage requirements as a muon can- 
didate track, and utilizes the same beam-constrained 
track fit. The track is required to have pr > 18 GeV, 
a kinematic region where the trigger track-finding ef- 
ficiency has no pt dependence. 

The clustering of showers in the CES produces an 
energy- weighted position at the electron shower max- 
imum. We require the CES cluster to be well sepa- 
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FIG. 14: The muon identification efficiency as a function FIG. 15: The electron track trigger efficiency as a function 
of the recoil component in the direction of the muon (itn). of track r\, for electrons identified in the calorimeter. The 

solid line shows the double-Gaussian parametrization of 

the data. 



rated from the edges of the towers, |CES x\ < 18 cm 
and |CES z\ > 9 cm. The cluster z position is com- 
pared to the extrapolated track z position, and the 
difference is required to be less than 5 cm, consis- 
tent with the trigger requirement. The ratio of the 
measured calorimeter energy to the track momentum, 
E/p, must be less than 2. 

Electrons are differentiated from hadrons by their 
high fraction of energy deposited in the electromag- 
netic calorimeter. The electron's EM energy is mea- 
sured in two neighboring towers in 77, while the energy 
collected in the hadronic calorimeter is measured in 
three towers. The ratio, E^ad) 'Eem , is required to 
be less than 0.1. Only the EM calorimeter measure- 
ment is used to determine the electron's px- 

An electron shower will typically be confined to 
a single tower, with a small amount of energy flow- 
ing into the nearest tower in 77. We define an error- 
weighted difference between the observed and ex- 
pected energies in the two towers neighboring the 
electron in the r\ direction [56| : 

L shr = 0.UJ2 , 1 . 1 (2°) 

i yjQ.U^E^ + {&Et xp y 

where E°f^ is the energy in a neighboring tower, E exv 
is the expected energy contribution to that tower, 
AE^ xp is the RMS of the expected energy, energies 
are measured in GeV, and the sum is over the two 
neighboring towers. We require Lshr ^ 0.3, consis- 
tent with the trigger criterion ( Section IIII A 4[) . 

The Z — > ee background is highly suppressed by 
the ut < 15 GeV requirement for the W boson sam- 



ple. Residual background results from electrons pass- 
ing through dead calorimeter regions, which reduces 
ut and increases p T . We remove events from the W 
sample if a track with pt > 20 GeV and \do\ < 0.3 
cm extrapolates to a calorimeter region with reduced 
response (|CES x\ > 22 cm or |CES z\ < 6 cm), 
and the track's calorimeter isolation is < 0.1 (Sec- 
tion IIV Ap . The full W — > ev selection results in a 
sample of 63,964 candidate events in (218.1 ± 12.6) 
pb _1 of integrated luminosity. 

The track selection in the single-electron trigger 
(Section [III A 4|) results in an ^-dependent trigger ef- 
ficiency for reconstructed electrons (Fig. [T5|) . We 
study this efficiency using W events selected with a 
trigger where the track requirements are replaced by 
a p T threshold. The efficiency decreases as \rj\ de- 
creases because the reduced path length reduces the 
ionization charge collected by each wire, thus reduc- 
ing the single hit efficiency. There is an additional 
decrease in efficiency due to the dead region at \z\ < 2 
mm. Electrons crossing this region at track |ry| = 
are not included in the efficiency plot, since we only 
measure electrons with |CES z\ > 9 cm. Thus, at 
1 7/ 1 = there is no inefficiency due to the dead COT 
region, and the measured efficiency increases. 

We measure the uu dependence of the electron 
identification efficiency (Fig. 1 16[) using Z — > ee 
events, selected with one electron passing the W elec- 
tron candidate criteria and a second "probe" elec- 
tron identified as an EM energy cluster with pt > 30 
GeV, an associated track with px > 18 GeV, and 
E/p < 2. Since the probe electron definition includes 
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FIG. 16: The electron identification efficiency measured 
in Z — > ee data as a function of the recoil component 
in the direction of the electron Background is sub- 

tracted using the number of like-charge lepton events ob- 
served at a given itn. The E/p < 2 requirement is not 
included in this efficiency measurement. 



lix, and q is the muon charge. The a priori mo- 
mentum scale is determined by the measurements of 
the magnetic field and the radius of the tracker. At 
CDF, eB/2 = 2.11593 x 10~ 3 GeV/cm, where B is 
measured using an NMR probe at a COT endplate. 
Measurements of the local field nonuniformities and 
tracker geometry were performed during construction 
and installation and are used to determine the posi- 
tions of individual track hits. We find these measure- 
ments provide an a priori momentum scale accuracy 
of « 0.15%. 

We refine the momentum scale calibration with 
data. Using reconstructed cosmic ray muon tracks, 
we align the relative positions of the tracker wires. 
Track-level corrections derived from W — > ev data 
reduce relative curvature bias between positive and 
negative particles. Finally, we perform an abso- 
lute calibration of the momentum scale using high- 
statistics data samples of J/ip, T, and Z boson de- 
cays to muons. The final calibration is applied as a 
relative momentum correction Ap/p to the W boson 
data and has an accuracy of w 0.02%. 



an E/p requirement, this cut is not included in the 
efficiency measurement. We instead study the unbi- 
ased E/p < 2 efficiency by recalculating E and it|| 
for towers separated in </> from the identified electron 
in W — > ev events, and find no significant uu de- 
pendence in this efficiency. In the simulation we use 
b = 0±0.54 x 10~ 3 , obtained by fitting the measured 
efficiencies to the function in Eqn. \T§\ 

We vary b by ±3<r in pseudoexperiments and as- 
sume linear variation of mw with b to derive uncer- 
tainties of Smw — 3, 5, and 16 MeV for the tot, Pt, 
and p T fits, respectively. Since b is measured with dif- 
ferent data samples for the electron and muon chan- 
nels, there is no correlation between the correspond- 
ing systematic uncertainties. 

V. TRACK MOMENTUM MEASUREMENT 

Muon momenta are determined from helical fits to 
tracks reconstructed using COT information. The 
momentum resolution of prompt muons is improved 
by constraining the helix to originate from the trans- 
verse beam position. A given muon's transverse mo- 
mentum is determined by the Lorentz equation, 

mv 2 / R — evB, ,_ , 

p T = eB/(2\c\), { } 

where B is the magnetic field, R is the radius of 
curvature, c = q/(2R) is the curvature of the he- 



A. COT Alignment 

The COT contains 30,240 sense wires for measur- 
ing the positions of charged particles passing through 
the detector. The position measurements rely on an 
accurate knowledge of the wire positions throughout 
the chamber. We determine these positions using a 
combination of alignment survey, computer model- 
ing, and cosmic-ray muon data. Any remaining bi- 
ases in track parameter measurements are studied 
with J /i\) — > fifi and W — > ev data, from which fi- 
nal track-level corrections are derived. 

After construction of the COT endplates, the po- 
sition of each 12- wire cell was measured with an ac- 
curacy of ±13 fim using a coordinate measuring ma- 
chine. The effect of the load of the wire plane and 
field sheets was modeled with a finite element analysis 
(FEA) and found to cause an endplate bend towards 
z = cm, with the maximum bend of ~ 6 mm in 
the fifth superlayer [29] . An equivalent load was ap- 
plied to the detector and further measurements found 
the FEA to be accurate to within w 20%. The FEA 
results were scaled to match the measurements, and 
the positions determined from the FEA were set as 
the directly-determined cell positions. 

While each cell position determines the average po- 
sitions of its 12 sense wires within the chamber, sev- 
eral effects create a non-linear wire shape as a func- 
tion of z. Gravity has the most significant effect, 
causing each wire to sag ~ 260 /im in y at z = 
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FIG. 17: The definitions of the local tilt (At) and az- 
imuthal shift (RAcj>) alignment corrections applied to 
each COT cell. 



cm. Electrostatic deflection towards the nearest field 
sheet occurs for cells where the sense wire is not cen- 
tered between the field sheets. By construction, the 
wires are slightly offset within a cell; in addition, the 
gravitational sag of the field sheets is larger than that 
of the sense wires, resulting in an electrostatic deflec- 
tion that partially counteracts the sag of the sense 
wires. Combined, the electrostatic effects cause a in- 
dependent wire shift that has a maximum of 74 /im 
at <j> = 145° and z = cm. The gravitational and 
electrostatic effects were combined to determine the 
best a priori estimate of the wire shapes. 

Starting from the predicted cell and wire positions, 
we develop in situ corrections based on cosmic-ray 
muon data taken during pp crossings with the sin- 
gle muon trigger. The data are selected by requiring 
exactly two reconstructed tracks in the event, elim- 
inating effects from overlapping hits from collision- 
induced particles. Since the two tracks on oppo- 
site sides of the COT result from a single cosmic-ray 
muon, we refit both tracks to a single helix and de- 
termine hit residuals with respect to this helix (57j . 
For each cell, we use the residuals to determine a 
tilt correction about its center, and a shift correction 
along the global azimuth (Fig. [T"7|) . We show the 
tilt and shift corrections for the inner superlayer of 
the west endplate in Fig. [HI after removing global 
corrections. We apply these corrections to each cell 
of each superlayer in each endplate. In addition, we 
measure a relative east-west shift and include it in 
each cell's correction. 

We combine the cell-based corrections with wire- 
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FIG. 18: The local tilt (top) and azimuthal shift (bottom) 
alignment corrections applied to each cell of the inner 
superlayer of the west endplate. Not shown are a global 
0.0021 tilt correction and a small global rotation and shift 
of the COT that does not affect track measurements. 



based corrections for the shapes of the wires between 
the endplates. We measure these corrections as func- 
tions of z and radius R using the differences in the 
measured do and curvature parameters for the helix 
fits on opposite sides of the COT for a cosmic ray 
muon. The corrections are applied as additional off- 
sets A£ of the wires at z = cm, with a parabolic 
wire shape as a function of z. The corrections include 
a radial dependence, 

A£ = -160 + 380(i?/140) - 380(i?/140) 2 , (22) 

where R is measured in cm and A£ in /im. Figure 
[T9l shows the gravitational and electrostatic shifts of 
a wire as a function of z at <j> = it, as well as the 
data-based correction at R = 130 cm (the outer su- 
perlayer) . 

The cell- and wire-based corrections are imple- 
mented for the track-finding and fitting stage, and re- 
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FIG. 19: The net wire shift in y as a function of z from 
gravitational sag only (solid), including electrostatic ef- 
fects (dashed), and including data-based corrections from 
Eq. ((22} (dotted). The shift is shown at (f> — T an d 
R = 130 cm. 



duce the measured hit resolution for high-momentum 
muons from wl80 fim to «140 fim,. Final track-based 
corrections are applied to the measured track cur- 
vature, which is inversely related to the transverse 
momentum [Eq. ([21])] . Expanding the measured cur- 
vature c as a function of the true curvature Ct in a 
Taylor series around zero, 



c = ei + (1 + e 2 )c t + £3Cf + e 4 cf + 



(23) 



the terms even in ct cause biases in positive tracks rel- 
ative to negative tracks, which tend to cancel when 
the two are averaged. The term linear in ct scales the 
true curvature and is determined by the momentum 
calibration. The e±c\ term is the first to directly af- 
fect mass measurements and is suppressed by the cf 
factor at low curvature (high momentum). 

Corrections for high-momentum tracks from W 
and Z decay particles are determined using the dif- 
ference in E/p for e + and e~ from W decays, which 
should be zero in the absence of misalignments. This 
difference can be used to constrain ex, the first term 
in the Taylor expansion. Figure [20] shows the differ- 
ences in E/p as functions of cot 9 and <j>, before and 
after corrections of the following form: 
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FIG. 20: The difference between e + and e~ E/p as a 
function of cot 6 (top) and <f> (bottom) before (triangles) 
and after (diamonds) track-level corrections. 



a relative rotation of the east and west endplates 
(aicotf?); and a mismeasurement of the beam posi- 
tion (bi s'm((j> + 0.1)). The measured values of the pa- 
rameters ao, ai, a2, b±, and 63, are shown in Table iLLTl 
Varying a\ by ±3<r in pseudoexperiments and as- 
suming linear variation of the momentum scale with 
ax, we find the ax uncertainty results in a relative mo- 
mentum scale uncertainty of ±0.07 x 10~ 3 for W and 
Z boson mass measurements. The other parameter 
uncertainties, as well as residual higher-order terms, 
have a negligible impact on the momentum scale for 
the mw measurement. 



6c = ao + ax cot 9 + a 2 cot 2 9+ 

bx sin(0 + 0.1) + 63 sin(30 + 0.5). 



(24) 



The terms can be interpreted as arising from the 
following physical effects: a relative rotation of the 
outer edge to the inner edge of each endplate (ao); 



B. J/lp — > jifi Calibration 

With a measured a x BR of 16.3±^;| nb [H, 
J/i/> mesons are the Tevatron's most prolific source 
of resonant decays to muon pairs. In addition to 
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Parameter Value (xlO 7 cm x ) 



a -0.66 ±0.17 
oi -1.6 ±0.3 

a 2 -2.1 ±0.5 

6i -2.1 ±0.2 

b 3 5.7 ±1.7 



TABLE III: The parameters used to correct the track 
curvature of electrons and muons from W and Z boson 
decays. The values and statistical uncertainties are deter- 
mined from fits to the E/p difference between positrons 
and electrons. 

its high statistics, the J/ijj's precisely known mass 
(m J/v , = 3096.88 ±0.04 MeV [58jn_and narrow width 
(r J/v , = 0.0934 ± 0.0021 MeV make it a key 
component of the track momentum calibration. We 
perform measurements of the J/ip mass as a function 
of mean inverse muon px to determine a momentum 
scale correction and extrapolate to the high-p-r region 
relevant for W and Z boson decays. 

1. Data Sample 

The J/ip data sample is collected with a Level 1 
trigger requiring one px > 1.5 GcV XFT track with a 
matching CMU track segment, and a second pt > 1.5 
(2) GeV XFT track with a matching CMU (CMX) 
segment. At Level 3, the two corresponding COT 
tracks must have opposite charge and consistent z 
vertex positions (|Azo| < 5 cm), and must form an 
invariant mass between 2.7 and 4 GeV. The resolution 
on the invariant mass measurement degrades at high 
track momentum, so to avoid trigger bias the mass 
range is extended to 2 GeV < m M/J < 5 GeV when 
the pt of the muon pair p T ^ is greater than 9 GeV. 

Candidate events are selected offline by requiring 
two COT tracks, each with p T > 2 GeV, |d | < 0.3 
cm, and > 7 hits on each of the eight superlayers. 
The tracks must originate from a common vertex 
(|Azo| < 3 cm) and form an invariant mass in the 
range (2.95, 3.21) GeV. 

A significant fraction (« 20%) of the J/ip mesons 
in our data sample result from decays of B hadrons, 
which have an average proper decay length of ~ 0.5 
mm. The muons from the J/ip decay can thus orig- 
inate outside the beam radius. Therefore, no beam 
constraint is applied in the COT track fit of muon 
candidates from J/ip decays. 

The total sample consists of 606,701 J/ip candi- 
dates in (194.1 ± 11.3) pb _1 of integrated luminosity. 



2. Monte Carlo Generation 

We use pythia [5!| to generate J/ip — ► /i^u events, 
from which templates are constructed to fit the data 
for the momentum scale. The shape of the m MM dis- 
tribution from J/ip decays is dominated by the pr- 
dependent detector resolution. We therefore model 
the distribution as well as the px and relative 
Pt of the muons in a J/ip decay. To obtain an ade- 
quate model, we empirically tune the generated J/ip 
kinematics to describe the relevant data distributions 
for the J/ip mass fits. 

To tune the p T ^ distribution, we boost the J/ip 
momentum by changing its rapidity (yj/^) along its 
direction of motion pj/^. In 50% of the generated 
events we multiply yj/^p by 1.215, and in the other 
50% we multiply it by 1.535. The decay angle 9* 
in the J/ip rest frame relative to pj/^ is tuned by 
multiplying cot 9* by 1.3. After tuning, the simula- 
tion matches the relevant background-corrected data 
distributions, as shown in Fig. 1211 

The PYTHIA event generator does not include en- 
ergy loss due to final-state photon radiation from the 
muons in J/ip decays. To simulate this effect, we 
scale each muon's momentum by a factor x deter- 
mined from the following leading-log probability dis- 
tribution for soft photon radiation [59 . l60| : 

f{x)=(3{l-xf-\ (25) 

with 

/3=^[ln(Q 2 K)-l] (26) 
and Q 2 = m 2 j^. 

3. Momentum Scale Measurement 

The momentum scale is calibrated using J/ip de- 
cays by fitting the dimuon mass as a function of mean 
inverse pt of the two muons, and then extrapolating 
to high pt {{pt 1 ) ~ GeV" 1 ). This procedure re- 
sults in a track momentum calibration accuracy of 
0.025%. 

The momentum scale calibration requires an accu- 
rate modeling of the muon ionization energy loss in 
the tracker. Each muon passing through the silicon 
and COT detectors loses on average 9 MeV at normal 
incidence. The combined effect on the reconstructed 
m w is about 0.6% of mj/^, a factor of » 20 larger 
than our total uncertainty. Since the ionization en- 
ergy loss Ej varies only logarithmically with pt (Sec- 
tion IIII B~Tjl . the relative effect on the reconstructed 
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FIG. 21: The J/ip —> rt/x data (points) and tuned simula- 
tion (histogram) distributions of pt^ (top) and 
(bottom). The J^^/Pt i s equal to the sum of the 
track curvatures of muons from a J/ip decay, divided by 
2.11593 x 10" 3 . 



mass is: 



2pg 



El 
2ptf 



Ej (p^). (27) 



Thus, in a linear fit of Am/m as a function of mean 
inverse px, a non-zero slope approximately corre- 
sponds to Ej. Since we model the ionization energy 
loss based on the known detector material, this slope 
should be zero. We however find that we need to 
scale down the ionization energy loss from the detec- 
tor parametrization (Section lHI B ip by 6% to achieve 
a zero slope. We show the result of this tuning in 
Fig. 1221 replacing Am/m on the y— axis with the rel- 
ative momentum correction Ap/p to be applied to the 



data in order to measure 



= 3096.88 MeV. The 



tuning is based on a (p? 1 ) region of (0.1, 0.5) GeV 1 , 
divided into eight bins. We find a scale correction of 
Ap/p = [-1.64±0.06(stat)] x 10~ 3 from a linear fit 
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FIG. 22: The fractional momentum correction for data as 
a function of the mean inverse momentum of the muons 
from J/ip decays. In a linear fit, the intercept corresponds 
to the scale correction relevant for W and Z boson decays, 
and the slope corresponds to the remaining unmodeled 
ionization energy loss after material tuning. The uncer- 
tainties are statistical only. 



to Ap/p as a function of (p^ 1 ). 

Each Ap/p value in Fig. [22] is extracted via a 
binned likelihood fit to the m MM distribution for each 
(p^ 1 ) bin. Since the mass resolution varies signifi- 
cantly with (p^ 1 ), the fit ranges are adjusted from 
3.08 ±0.13 GeV for (p^) = (0.1,0.15) GeV" 1 to 
3.08 ±0.08 GeV for (p^ 1 ) = (0.45,0.5) GeV" 1 . The 
background is modeled as a linear function of m^, 
with normalization and slope determined from upper 
and lower sideband regions whose combined width 
is equal to that of the mass fit window. The re- 
sults of the fits in the (p T l ) = (0.15,0.2) GeV" 1 
and (p T l ) = (0.25,0.3) GeV" 1 ranges are shown in 

Fig. m 

The J/ip momentum calibration includes correc- 
tions to the curvature c derived from the measured 
dimuon mass as a function of A cot 9 between the 
positive and negative muons from the J/ip decay. Bi- 
ases linear in A cot 9 are removed with a curvature 
correction linear in cot 9: 



Sc= [(-7±1) x 10 



" 7 cm" 1 ' 



cot I 



(28) 



where the uncertainty is statistical only. Biases 
quadratic in A cot 9 are removed with the following 
correction to the absolute length scale of the COT 
along the z axis (statistical uncertainty only): 



5 cot 9 = [(-3.75 ± 1.00) x 10" 4 ] coti 



(29) 
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FIG. 23: The m^/j, fits to data (points) with (py 1 ) = 
(0.15,0.2) GeV- 1 (top) and (p" 1 ) = (0.25,0.3) GeV' 1 
(bottom). The best fit to the m MM = (3.08 ± 0.12) 
GeV (top) and m MM = (3.08 ± 0.11) GeV (bottom) 
regions correspond to momentum scale corrections of 
(-1.54 ± 0.09) x 10~ 3 (top) and (-1.65 ± 0.04) x 10" 3 
(bottom). The arrows indicate the fit regions and the 
uncertainties are statistical only. The fit \ 2 can be im- 
proved by adjusting the final-state radiation model, and 
this effect is incorporated into the systematic uncertainty 
(Section |VB4|. 



The cot 61-dependent correction to the curvature 
[Eq. (|2"8))] is larger than the correction derived from 
E/p in W — > ev data (oi in Table HIT)). Muons from 
J / ip decay have a broader curvature range and thus a 
greater dependence on misalignments affecting higher 
order terms in curvature. Since we derive a curvature 
correction averaged over all of the terms in Eq. (|2"5|) , 
the J/ip correction can be larger than the correction 
for electrons and muons from W and Z boson decays. 



Systematic uncertainties on the momentum scale 
correction extracted from J/t/j — ► /i/i decays (Ta- 
ble IIV[) are dominated by the incompleteness of the 
QED and energy loss models. At low muon px (high 
(pt 1 ))' the mass fits become increasingly sensitive to 
QED and energy loss modeling because of the better 
mass resolution and higher statistics. Since we only 
model the mean ionization energy loss, our modeling 
of the mass region below the peak is imperfect. Addi- 
tionally, our neglect of higher-order QED corrections 
affects the modeling of this region. We study possi- 
ble bias from our incomplete model by changing the 
Q 2 value in the photon radiation probability function 
[Eq. I|25[)] such that the x 2 of the inclusive m PM fit is 
minimized. We find that this change affects Ap/p by 
0.2 x 1CT 3 . 

If there is a relative tilt between the solenoid and 
the tracker axes, the extracted momentum scale cor- 
rection will have a linear dependence on A cot 9. In 
addition, incomplete corrections of the magnetic field 
nonuniformities near the ends of the solenoid can 
cause a quadratic cot 9 variation. We study the 
cott9 dependence of Ap/p using J ftp decays where 
both muons are measured in the same cot 9 region 
(| A cot 6 (nfJ>)\ < 0.1). We find that if we correct 
for the observed quadratic dependence, the extracted 
Ap/p changes by 0.1 x 10~ 3 . 

The uncertainty on the material correction propa- 
gates to a momentum scale uncertainty of 0.06 x 10 -3 
when extrapolated to high momentum, as shown 
in Fig. [2H An additional statistical uncertainty of 
0.01 x 10~ 3 on the scale is determined by fixing the 
material correction and fitting for the scale. 

The statistical uncertainties on the J ftp alignment 
corrections [Eq. and ([29])] have a 0.05 x 10" 3 



effect on Ap/p. We test our model of the m MM line- 
shape by changing the fit range by ±20%, and find a 
±0.05 x 10~ 3 change in Ap/p. 

We apply the same pr thresholds offline as in the 
trigger for muons with CMU segments. Since we do 
not model a j?T-dependent trigger efficiency, any inef- 
ficiency could cause a bias in the reconstructed m Mfl . 
We investigate this possibility by varying the offline 
Pt thresholds by ±5%, and find a Ap/p variation of 
±0.04 x 10~ 3 . 

The quality of the fit is highly sensitive to the hit 
resolution model, but the momentum scale correction 
is not. Changing the simulated COT hit resolution 
by ±10 /Ltm, which corresponds to a > 10ct statis- 
tical variation, results in a ±0.03 x 10 -3 change in 
Ap/p. We include this in our systematic uncertainty 
estimate. 
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A ±0.03 x 10~ 3 uncertainty on Ap/p from the back- 
ground model is determined by changing its linear de- 
pendence on to w to a constant. Finally, the world- 
average J/tp mass value used in this measurement 
contributes ±0.01 x 10~ 3 to the uncertainty on Ap/p. 

The final momentum scale correction derived from 
J/ijj data is: 

Ap/p = (-1.64 ± 0.25) x 10" 3 . (30) 



C. T — > n^i Calibration 

The bb resonance T provides a complementary mo- 
mentum scale calibration tool to the J/tp. Its pre- 
cisely measured mass mx = (9460.30 ± 0.26) MeV 
[rl| is three times larger than that of the J/tp, so 
an T momentum scale calibration is less sensitive to 
the material and energy loss model than that of the 
J/tp. Because the bb resonances are the highest mass 
mesons, long-lived hadrons do not decay to the T and 
the muons from T decay effectively originate from 
the collision point. We improve the accuracy of the 
muon measurements by constraining their tracks to 
the beam position, which is the same procedure ap- 
plied to the W and Z decay lepton tracks. 

The T data sample is based on the same Level 
1 trigger as the J/ip sample (Section IVB The 
Level 3 requirements are: one reconstructed track 
with p T > 4 GeV and matching CMU and CMP 
track segments (CMUP); a second track with oppo- 
site charge to the first, px > 3 GeV, and a matching 
CMU or CMX track segment; and a reconstructed 
mass of the two tracks between 8 and 12 GeV. Of- 
fline, the pt thresholds are increased to 4.2 (3.2) GeV 
for the track with a CMUP (CMU or CMX) track seg- 
ment, and each track must have \do\ < 0.3 cm and at 
least 5 hits in at least 3 axial and 3 stereo supcrlay- 
ers. The two tracks are required to have a common 
vertex (|Az | < 3 cm). 

We model T production and decay using PYTHIA 
[5^|, to which we apply the same tuning procedure 
as for J/tp generation. The data pj, distribution is 
matched in simulation by boosting the rapidity of 
each decay muon by 0.07?/t along pxj where yr is the 
T rapidity. Radiation of photons from the final state 
muons is simulated using the probability distribution 
of Eqns. [25l and [26l The pt^ distribution is shown in 
Fig. [2H after subtracting background from the data. 

We test any possible beam-constraint bias by sep- 
arately reconstructing charged muon tracks from T 
decays with and without incorporating the beam 
constraint. For the sample with beam-constrained 
tracks we fit for tot in the region 9.28 GeV < 
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FIG. 24: The p 1 ^ 1 (top) and p^. (bottom) distributions 
from T — > fifi decays for the data (points) and simulation 
(histogram) . 



m uu < 9.58 GeV, while for the sample with non- 
beam-constrained tracks we fit the region 9.25 GeV < 
to w < 9.61 GeV. In (190.8 ± 11.1) pb" 1 of inte- 
grated luminosity, we have 34,618 T candidates with 
beam-constrained tracks and 35,622 candidates with 
non-beam-constrained tracks. The two momentum 
scale measurements are shown in Fig. [25] and are con- 
sistent at the 2a level when correlations are taken 
into account. We define the T result to be the 
mean of the two values, and take half their difference 
(Ap/p — 0.06 x 10 -3 ) as a systematic uncertainty on 
the measurement. 

The remaining systematic uncertainties on the mo- 
mentum scale measurement with T decays are com- 
mon to those of the measurement with J/ip decays. 
We use the same procedures as with the J/tp calibra- 
tion to estimate the sizes of the uncertainties, with 
one exception. Since the T sample has < 10% of the 
statistics of the J/ip sample, the QED and energy loss 
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Source 


J/ip (xlO~ 


" 3 ) T (xlO -3 ) 


Common (xl0~ 3 ) 


QED and energy loss model 


0.20 


0.13 


0.13 


Magnetic field nonuniformities 


0.10 


0.12 


0.10 


Beam constraint bias 


N/A 


0.06 





Ionizing material scale 


0.06 


0.03 


0.03 


COT alignment corrections 


0.05 


0.03 


0.03 


Fit range 


0.05 


0.02 


0.02 


Trigger efficiency 


0.04 


0.02 


0.02 


Resolution model 


0.03 


0.03 


0.03 


Background model 


0.03 


0.02 


0.02 


World-average mass value 


0.01 


0.03 





Statistical 


0.01 


0.06 





Total 


0.25 


0.21 


0.17 



TABLE IV: Uncertainties on the momentum scale correction derived from the J/ip and T mass measurements. 



model cannot be tested with the \ 2 °f the T — > /i/i 
mass fit. Instead, we change Q in the photon radia- 
tion probability by the amount estimated for the J/ip 
systematic uncertainty ( Section I V B 4[) . We find that 
this variation affects Ap/p by ±0.13 x 10~ 3 in the T 
calibration. 

The final result of the T calibration is: 

Ap/p = (-1.44 ± 0.21) x 10~ 3 . (31) 

We have verified that this result has no time de- 
pendence, at the level of the statistical precision of 
±0.13 x 10~ 3 . When combined with the momentum 
scale correction from the J /if) calibration, we obtain: 

Ap/p = (-1.50 ± 0.19) x 10~ 3 . (32) 

D. Z — > /i/i Calibration 

Given the precise momentum scale calibration from 
the J/?/ 1 and T decays, we measure the Z boson mass 
and compare it to the world-average value mz = 
(91187.6 ± 2.1) MeV [ll|. We then use the world- 
average mz to derive an additional Ap/p calibration 
and combine it with that of the J/ip and T decays. 

The systematic uncertainties of the mz measure- 
ment are correlated with those of the mw mea- 
surement, so a momentum scale calibration with Z 
bosons can reduce systematic uncertainties on the 
mw measurement. However, the statistical uncer- 
tainty from the Z — > fi[i sample is significantly larger 
than the calibration uncertainty from J/ip and T de- 
cays. Thus, the main purpose of the mz measure- 
ment is to confirm the momentum scale calibration 
and test our systematic uncertainty estimates. 



The Z boson data sample is selected using the same 
single-muon trigger and offline muon selection as for 
the W boson sample (Sections !"!!! A 4l and lIV A|) . with 
the exception that we remove the requirement of a 
track segment in a muon detector for one of the 
muons from the Z boson decay. Removing this re- 
quirement significantly increases detector acceptance 
while negligibly affecting background. Z boson can- 
didates are defined by 66 GeV < m^ < 116 GeV, 
Pjt 1 < 30 GeV, \At (fi,fi)\ < 3 ns, and oppositely 
charged muons. A muon track's to is defined as the 
time between the pp bunch crossing and the muon's 
production, and should be (0 ± 1) ns for Z — > fifi 
production and decay. The track to is measured us- 
ing the time information from the track hits in the 
COT by incorporating to into the helical fit. The 
|Afo| < 3 ns requirement effectively removes cosmic 
ray muons passing through the detector. An addi- 
tional cosmic ray identification algorithm reduces 
this background to a negligible size. After applying 
all selection criteria, the Z — > sample contains 
4,960 events in (190.8 ± 11.1) pb" 1 of integrated lu- 
minosity. 

We model Z boson production and decay using the 
RESBOS [6l| event generator and a next-to-leading 
order QED calculation of photon radiation from the 
final-state muons [13] (Section IIX|) . For m w near 
the Z boson resonance, the photon propagator and 
Z/7* interference make small contributions to the 
shape of the m w distribution. We separately sim- 
ulate these components and include them as fixed 
"background" to the Z lineshape. We measure mz 
using a binned likelihood template fit to the data in 
the range 83 GeV < m w < 99 GeV (Fig. [2SJ). Our 
measurement of m z = [91.184 ± 0.043(stat)] GeV is 
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FIG. 25: The momentum scale correction Ap/p derived 
from binned likelihood fits to the m MM data distribution 
(points) in the region dominated by T — > /i/i decays. The 
small difference between fits using tracks with (top) and 
without (bottom) a beam constraint is incorporated into 
the systematic uncertainty. The arrows indicate the fit 
region and the uncertainties are statistical only. 



in good agreement with the world-average value of 
m z = (91.188 ± 0.002) GeV 

Systematic uncertainties on mz are due to the mo- 
mentum scale calibration (17 MeV), alignment cor- 
rection uncertainties (7 MeV), and incomplete mod- 
eling of higher-order QED corrections (14 MeV). The 
combined statistical and systematic uncertainty is 49 
MeV. 

Given the precise world-average measurement of 
mz, we use the Z boson resonance as an additional 
calibration input to Ap/p. We find that adding the 
mz information reduces Ap/p and its uncertainty by 
less than 0.01 x 10 -3 each. 

Incorporating the alignment uncertainty (Sec- 
tion IV Aj) into Ap/p from Eqn. [32] gives the mo- 
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FIG. 26: The m z fit to data (points) in the 83 GeV 
< rrifi^ < 99 GeV region (denoted by arrows). The un- 
certainty is statistical only. 



mentum scale correction applicable to the W boson 
sample: 

&p/p = (-1.50 ±0.21) x 10" 3 . (33) 

The corresponding uncertainty on the mw fits in the 
muon channel is 17 MeV. 



VI. ELECTRON ENERGY MEASUREMENT 

An electron's energy is measured from its shower 
in the electromagnetic calorimeter. We perform an 
initial data calibration by scaling the measured en- 
ergy such that a Gaussian fit to the reconstructed 
dielectron mass in a region dominated by Z decays 
(86-98 GeV) gives a mean of 91 GeV. This is slightly 
below the world-average mz because the Gaussian fit 
is biased by the energy lost to final-state photon ra- 
diation (Section IIX D[) . This initial data calibration 
is accurate to sa 0.15%. 

To model the data, the simulated calorimeter en- 
ergy is scaled to match the measured E/p distribu- 
tion of electrons in W — > ev events. A calibrated 
data measurement would result in an E/p of unity 
for electrons that do not radiate before entering the 
calorimeter, and deposit all of their energy in the 
EM calorimeter. We verify that the E/p calibration 
is unbiased by using it to measure mz in dielectron 
events. Given consistency of the measured mz with 
the world-average value, we incorporate the mz fit 
into the calibration. The final calibration has an ac- 
curacy of 0.037%. 
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A. E/p Calibration 

We transfer the precise tracker calibration to the 
calorimeter using the ratio of electron calorimeter 
energy to track momentum, E/p. The material 
from the beam pipe to the inner COT wall causes 
bremsstrahlung that affects the measured position of 
the E/p peak, and this material is scaled such that 
the simulation matches the data in the high E/p re- 
gion. The non-linearity of the energy scale is removed 
by applying a correction to the simulation scale as a 
function of the calorimeter shower pt [Eq. (|17p]. Fi- 
nally, corrections are applied to the data to improve 
uniformity in response as a function of detector tower 
and time. After the complete set of corrections and 
simulation calibrations, the simulation energy scale 
Se is determined from a maximum likelihood tem- 
plate fit to the E/p peak region. 

The shape of the E/p distribution has a strong 
dependence on the material upstream of the COT. 
Bremsstrahlung in this material reduces the mea- 
sured electron momentum in the tracker while leaving 
the measured calorimeter energy unchanged, since 
photons are radiated collinearly with the electron and 
deposit their energy in the same calorimeter tower as 
the electron. Thus, the effect of bremsstrahlung is to 
shift the measured E/p to values > 1. If the material 
were not well modeled, the energy scale calibration 
would be biased to compensate for the mismodcling. 

A detailed accounting of the silicon and COT 
tracker material was performed at installation. In 
the early data-taking period, the radial distribution 
of photon conversions was compared between data 
and a full GEANT simulation. The amount of copper 
cable was increased by a few percent of Xq in the 
GEANT simulation to correct observed discrepancies, 
and the three-dimensional lookup table of material 
properties (Section IIII B 1[) was produced from this 
corrected GEANT simulation. 

For a final material tuning, we compare our 
parametrized simulation to the data in the high E/p 
region (1.19 < E/p < 1.85) of electrons from W bo- 
son decays. Using the region 0.85 < E/p < 1.19 for 
normalization, we perform a maximum likelihood fit 
to the 1.19 < E/p < 1.85 region in two bins (Fig. [27]) 
and measure a radiation length multiplicative correc- 
tion factor of S mat = 1.004 ± 0.009(stat) [63]. As a 
further consistency check of the material lookup ta- 
ble, we determine S ma t as a function of tower |^|, and 
find no statistically significant dependence on \r]\. 

Our simulation of electron interactions in the 
tracker and calorimeter accounts for most of the en- 
ergy dependence of the energy scale. Any residual 
non-linearity is incorporated as a per-particle cor- 



' 1 ' 5 E/p(W^ev) 

FIG. 27: The maximum likelihood fit to the tracker radi- 
ation length correction factor using the high E/p region 
(1.19 < E/p < 1.85) in W -> ev data (points). 



rection in the simulation (Section IIII B 2p . To mea- 
sure this non-linearity, we fit the E/p peak region 
(0.93 < E/p < 1.11) for the energy scale in bins of 
measured electron calorimeter pt (Fig. [55]) ■ The re- 
sulting energy scale measurements are fit as a linear 
function of pr, fixing the scale to 1 at the W boson 
sample's (pf^) = 39 GeV. The error-weighted aver- 
age, £ = [6± 7(stat)] x 10~ 5 , of the measurements of 
the non-linearity parameter from the W and Z boson 
samples is used in Eq. ([17]) . The linear fits in Fig. |2"5I 
where the simulation includes this correction, show a 
constant energy scale [64j . 

To improve the energy resolution of the data, we 
apply time-dependent and tower-dependent calibra- 
tions derived from low-energy EM clusters. At level 3 
the relevant trigger requires calorimeter and track pt 
greater than 8 GeV each, as well as electron iden- 
tification based on track-calorimeter matching and 
calorimeter shower shape properties. Offline, can- 
didates are required to have Had/EM< 0.05 and 
E + p > 22 GeV to remove any trigger bias. Us- 
ing the mean of the E/p range 0.8 — 1.25, we apply 
relative corrections of 0(3%) to remove variations as 
functions of tower and time. 

Because of bremsstrahlung radiation in the tracker, 
the mean E/p correction has a small bias that de- 
pends on the electron path length. Since the path 
length increases as |?7| increases, we perform a final 
^-dependent calibration of the data. Using template 
fits to the E/p peak region of the W — > ev sample in 
bins of | ?y |, we derive a relative correction for each bin. 
This calibration removes ~ 2% residual variation in 
the calorimeter energy response. 
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FIG. 28: The energy scale as a function of measured elec- 
tron calorimeter pr for W (top) and Z (bottom) boson 
decays. The fast simulation incorporates a per-particle 
non-linear response correction [Eq. (j!7p ] . The combina- 
tion of the linear fits to £ results in no energy dependence 
of the energy scale. 



FIG. 29: The energy scale calibration using the peak E/p 
region (0.93 < E/p < 1.11, denoted by arrows) in W — > 
ev data (points). 
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With the complete set of corrections applied to the 
data and simulation, we calibrate the simulation en- 
ergy scale using W — > ev events. The fit for Se 
[Eq. ([P?])] to the E/p peak region (Fig. [29]) has a sta- 
tistical uncertainty of 0.025%. Including systematic 
uncertainties due to S ma t (±0.011%) and the tracker 
momentum scale (±0.021%), we obtain a total uncer- 
tainty of 0.034% on the E/p calibration of the elec- 
tron energy scale. 

The E/p calibration requires an accurate simula- 
tion of electron radiation in the tracker. We test the 
track simulation by measuring mz (Section lVI Bp us- 
ing electron track information only. The measure- 
ment is a binned likelihood fit to the region 75 GeV < 
m ee < 99 GeV (Fig. |3"0|) , with mz as the fit param- 
eter. Because of the significant radiated energy loss, 



FIG. 30: The m z fit to data (points) in the 75 GeV 
< rn.ee < 99 GeV region (denoted by arrows), using re- 
constructed track information only. The uncertainty is 
statistical only. 



the test is less precise than the measurement using 
the calorimeter (Fig. l3"Tj) . Nevertheless, we obtain 
good consistency with the world-average mz , verify- 
ing that we do not have any significant mismodeling 
of electron radiation in the tracker. 



B. Z — > ee Calibration 

Using the _E/p-based calorimeter energy calibra- 
tion, we measure the Z boson mass from its decay 
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FIG. 31: The m z fit to data (points) in the 81 GeV 
< m ee < 101 GeV region (denoted by arrows). The un- 
certainty is statistical only. 



to two electrons. After confirming consistency of the 
result with the world-average mass, we fix mz to this 
value and produce a combined calibration from the 
electron E / p-based method and Z — > ee mass mea- 
surement. 

We select Z bosons using the same single-electron 
trigger and offline electron selection as for the W 
boson sample (Sections IIII A 41 and IIV B[) . and de- 
fine candidates as oppositely-charged electrons with 
66 GeV < m ee < 116 GeV and pf < 30 GeV. The Z 
boson sample contains 2,919 events in (218.1 ± 12.4) 
pb _1 of data. 

The sample includes a small component of multijet 
and W+ jet background. From a comparison of the 
data with like-sign electrons to a prediction of the 
full GEANT simulation, we estimate the background 
fraction to be < 0.5%. Since (m ee ) of the background 
is w 2 GeV less than that of the Z boson sample in 
the fit region, we estimate any corresponding bias on 
the measured mz to be < 10 MeV. 

The model for Z boson production and decay to 
electrons is the same as for the muon decay channel 
(Section IVDp . We use the resbos [6l| event gen- 
erator and a next-to-leading order QED calculation 
of photon radiation from the final-state electrons [62| 
( Section HX|) . We include the virtual photon exchange 
and Z/j* interference contributions as fixed "back- 
grounds" to the Z boson lineshape, and determine 
mz from a binned likelihood fit to the data in the 
range 81 GeV < m ee < 101 GeV (Fig. [51]). 

Systematic uncertainties on the mz measurement 
result from the E/p calibration (29 MeV), calorime- 
ter non-linearity measurement (23 MeV), and higher- 



order QED radiation (14 MeV). The measured mz — 
[91.190 ± 0.067(stat)] GeV is consistent with the 
world-average value m z = (91.188 ±0.002) GeV [HI], 
given the total uncertainty of 78 MeV on the mea- 
surement. 

The uncorrelated uncertainties in the combination 
of the mz and E/p calibrations are the uncertainty 
on the non-linearity parameter £, the statistical un- 
certainty on the mz measurement (0.073%), and the 
uncertainty on the E/p calibration (0.034%). Since 
the myy fit relies predominantly on the shape of the 
Jacobian edge of the distribution, the relevant 
electron transverse energies are in the « 40 — 45 GeV 
range. The uncertainty on the energy dependence of 
the scale from the Z boson mass is negligible, as the 
(pj) is about 42 GeV in this sample. The E/p-haseA 
calibration involves an extrapolation from {p c T ) = 39 
GeV, so it receives an additional uncertainty contri- 
bution of 23 MeV to the mw measurement from the 
non-linearity parameter £. Combining the two cali- 
brations, we obtain a total electron energy measure- 
ment uncertainty of 30 MeV on the mw measurement 
in the electron channel. Of this uncertainty, we take 
17 MeV to be 100% correlated with the muon channel 
through the momentum scale uncertainty. 



VII. RECOIL MEASUREMENT 

The recoil ut (Fig. [5]) in a W boson event results 
from quark or gluon radiation in the initial state, and 
from photon radiation in the initial and final states. 
A quark or gluon typically fragments into multiple 
hadrons, which are detected in the calorimeter. Ad- 
ditional energy from the underlying event is also mea- 
sured in the calorimeter and obscures the recoil mea- 
surement. Rather than rely on detailed modeling of 
the underlying event, we develop an empirical model 
of the recoil ut using Z boson events, where the four- 
momentum of the Z boson is measured precisely us- 
ing its leptonic decays. The model of the recoil energy 
measurement is tuned with these decays and applied 
to W boson events. 

We measure the recoil energy using all calorimeter 
towers except those with ionization or shower energy 
from the charged leptons. To reduce potential bias 
and facilitate our model parametrization, we correct 
the measured energy in each tower for acceptance dif- 
ferences resulting from an uncentered beam. In addi- 
tion, we improve the measurement resolution by cor- 
recting for response differences between the central 
and plug calorimeters. 
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FIG. 32: The maximal energy correction factors applied 
to each tower as a function of 77 j. The decreasing curve 
corresponds to towers closest to the beam line, and the 
increasing curve corresponds to towers on the opposite 
side. The correction factors reduce the cf> variation in 
acceptance due to the beam line displacement from the 
detector axis. 



A. Data Corrections 

The data used in this analysis have a relative offset 
of about 4 mm between the beam line and the center 
of the CDF II detector. This results in a variation in 
calorimeter acceptance as a function of <fr such that 
the calorimeter towers closest to the beam line have 
a larger acceptance, and thus a larger average mea- 
sured energy per tower. The variation is largest in the 
forward region, where the towers are in closest prox- 
imity to the beam line. We suppress this azimuthal 
energy variation by applying a threshold on the com- 
bined EM and hadronic tower pt of 5 GeV for tow- 
ers with detector \rj\ > 2.6. The threshold strongly 
suppresses the forward tower energy variation in W 
boson events, while retaining the energy from high- 
er hadronic jets in multijet events. All other towers 
have EM and hadronic energy thresholds of 20 MeV 
each. 

We reduce the residual azimuthal energy variation 
by applying a multiplicative correction factor to each 
measured tower energy according to the following em- 
pirical function (Fig. [32)) : 

Staler = 1 - . 6 (0 . 32 1 77 1 ) 4 • 74 sin(0 - 0.47). (34) 

This correction is determined using events collected 
by a minimum bias trigger, which requires evidence 
of an inelastic pp collision (Section IIII A 5|) . 

The relative energy scale between the central and 



forward calorimeters is initially determined from the 
calibration of high-p^ hadronic jets. The relative re- 
sponse has a significant energy dependence, however, 
and the initial calibration is not optimized for the low 
Pt particles relevant to the W boson recoil measure- 
ment. Using the E/p distribution of charged pions 
from minimum bias events, we find that a relative 
energy scale of « 12% between central and forward 
calorimeters is appropriate for particles with px < 2 
GeV, the momentum region of a typical recoil par- 
ticle. To maintain the mean recoil energy scale, we 
scale the central (forward) calorimeter tower energies 
up (down) by 5% (7%). This calibration improves the 
recoil resolution, and thus the statistical precision of 
the mw fits. It also minimizes the sensitivity of the 
recoil model to differences in phase space sampled by 
the selected W and Z boson decays. 



B. Lepton Tower Removal 

The recoil ut is measured as the sum of corrected 
Pt in all calorimeter towers (Sec. IVII Al) . excluding 
the towers in which the lepton(s) deposit energy. The 
exclusion of these towers also removes some recoil en- 
ergy from the measurement, thus causing a bias in u\ i . 
We measure this bias from the data and incorporate 
it in the simulation. 

An electron shower typically distributes energy to 
two calorimeter towers, but can also contribute to a 
third tower if the electron is near a tower edge. We 
remove each tower neighboring the electron's tower, 
as well as the corner towers closest to the electron's 
CES position (Fig. [33]). A muon near a tower edge 
can cross two towers, so we remove the two towers in 
-q neighboring the muon's tower (Fig. |3"4"|) . The tower 
window definitions are motivated by the presence of 
excess energy in a given tower above the background 
energy from the underlying event. 

We estimate the recoil energy flow into the ex- 
cluded towers, denoted by Aiti , using equivalent win- 
dows separated in <fi from the lepton in W — > Iv 
events. When simulating a W or Z boson event, we 
correct the simulated ut by a Ait| taken from the 
measured distribution. The simulated Au|| incorpo- 
rates its measured dependence on u\\ and u±, and lep- 
ton \rj\. These dependencies are shown for W — > [iv 
events in Fig. [35] and similar functions are defined 
for electrons. The incorporation of these functions 
preserves (Am||), which is 269 MeV for electrons and 
112 MeV for muons (with negligible statistical uncer- 
tainty) . 

To estimate the systematic uncertainty associated 
with modeling the tower removal, we study the varia- 
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FIG. 33: The average energy collected in the electromag- 
netic (top) and hadronic (bottom) calorimeters in the 
vicinity of the electron shower in W boson decays. The 
differences Acj> and A 7/ are signed such that positive differ- 
ences correspond to towers closest to the electron position 
at shower maximum. The central seven towers inside the 
box are removed from the recoil measurement. Statistical 
uncertainties on the values outside the box are 0(1 MeV). 
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FIG. 34: The average energy collected in the electromag- 
netic (top) and hadronic (bottom) calorimeters in the 
vicinity of the muon in W boson decays. The differences 
A<f> and A77 are signed such that positive differences cor- 
respond to towers closest to the muon position at the 
shower maximum detector. The central three towers in- 
side the box are removed from the recoil measurement. 
Statistical uncertainties on the values outside the box are 
0(1 MeV). 



tion of Am 1 1 in the data as a function of the (f> separa- 
tion from the lepton of the equivalent tower window. 
We take half the variation as a systematic uncer- 
tainty: 8 (5) MeV for removed electron (muon) tow- 
ers. To confirm our estimate of this uncertainty, we 
remove an additional window azimuthally opposite 
to the lepton (A\<fi = n), incorporate its model into 
the simulation, and compare the resulting simulation 
and data M|| distributions. We find the differences to 
be consistent within our quoted uncertainties. 



C. Recoil Model Parametrization 

The recoil consists of three separate components: 
radiation in the W or Z boson production; radiation 
from the spectator partons; and energy from addi- 
tional pp collisions in a given bunch crossing. We use 
the RESBOS [6l[ generator to predict the net pr distri- 



bution of radiation in the W or Z boson production, 
and minimum bias data for the px distribution from 
spectator partons and additional interactions. The 
parameters for the detector response to the recoil are 
measured in Z boson events. 

To facilitate tuning of the recoil model, we define 
axes such that quark and gluon radiation lies pre- 
dominantly along one axis, denoted as the "77" axis 
(Fig. [35]). This axis is chosen to be the angular bi- 
sector of the two leptons, whose angles are precisely 
measured. The orthogonal axis is denoted as the 
axis. 



1. Recoil Energy Scale 

We tune the simulation to match the observed de- 
tector response to the recoil radiation. The recoil re- 
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FIG. 35: The scales applied to the energy in the three 
removed muon towers in the simulation (solid lines), as 
functions of uu (top), u±_ (middle), and rj (bottom). The 
points show the scales measured using towers separated in 
<f> from the muon in W — * jw data. The scaling functions 
for removed electron towers have similar shapes. 
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FIG. 36: The definitions of the rj and £ axes in Z boson 
events. The quark and gluon radiation from the boson 
production points predominantly in the — ij direction. 



stants determined from the data. Figure [37| shows 
—Urf neas ■fij!'* 1 /Pt^i which approximates R, for Z bo- 
son decays to muons. The response R is less than 1 
due to calorimeter energy loss from particles curling 
in the tracker, particles passing through calorimeter 
cracks, and non-linearity of the hadronic calorimeter 
response. 

Projecting the lepton momenta and the recoil along 
the rj axis to obtain p 1 ^ and u v , the sum pz+u^ is sen- 
sitive to R. This sum is zero for R = 1, and positive 
for R < 1. We measure A = 0.635 ± 0.007(stat) and 
B = 6.68 ± 1.04(stat) by minimizing the combined 
X 2 of the electron and muon (p** + u n ) distributions 
as a function of (Fig. l38|). We determine A and B 
with the (pz +u v ) distribution rather than the distri- 
bution of Fig. [37] because (p" + u v ) is well-defined as 
— > GeV, while R is not. The parameters A and 
B are statistically uncorrelated by construction. We 
apply R(A, B) to the generated recoil ut in simulated 
W and Z boson events. 



2. Spectator and Additional pp Interactions 



sponse is defined as R = u™ eas ■ u^ rue /u^f ue , where 
Ur£ rue — —pj? is the generated net pr of the initial- 
state-radiation, and Urp s is the reconstructed vec- 
tor of this transverse momentum. 

To simulate the measured recoil, we parametrize 
the response as 

R(A,B) = A\n(ut£ ue + B)/ \n(15 + B), (35) 

where u*rf ue is in units of GeV, and A and B are con- 



The net pr fl° w from spectator quarks and ad- 
ditional interactions is negligible due to momentum 
conservation. However, detector resolution causes 
its measurement to generally be non-zero. The res- 
olution is predominantly determined by the energy 
sampling in the calorimeter, and we expect it to in- 
crease as the square root of the scalar sum J2pt of 
the calorimeter tower px- We plot the width of the 
fl T distribution, projected along the x and y axes, as 
a function of the ^2pt in minimum bias data. We 
parametrize the dependence as a power law, with the 
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FIG. 37: The mean -u T meas ■pf/p!]?', which approxi- 
mates the recoil response 7?, as a function of dimuon pr 
for m^fj, in the Z boson mass region. The simulation 
(solid) uses parameters fit from the electron and muon 
p" + Uri distributions, and models the data (circles) well. 



fitted result: 

x 0.5333 

a x , y = 0.3842 (J2pt) GeV, (36) 

where ^2pt is defined in units of GeV. The distri- 
bution of ^2pt from additional interactions, denoted 
Pmb, is parametrized as (Fig. |3"!)|) : 

Pub (j» oc i^p T )^ e" ^ 19 ' 98 , (37) 

with constants obtained from a fit to the minimum 
bias data. In our simulation, we draw a value of 
from this distribution, for the fraction of events 
containing at least one pp collision beyond that pro- 
ducing the W or Z boson. This fraction is calcu- 
lated from the average instantaneous luminosity of 
2.137 (2.014) x 10 31 cm" 2 s" 1 for W and Z boson 
data in the muon (electron) channel, and the as- 
sumed instantaneous luminosity per additional col- 
lision (3.3 x 10 31 cm~ 2 s" 1 ). 

The observed ^2pt from spectator partons in the 
pp — > W or Z boson interaction is modeled from the 
minimum bias data, which correspond to one or more 
pp collisions. We deconvolute the J^Pt spectrum of 
Eqn. [37J with the distribution of the number of col- 
lisions in minimum bias data to derive the following 
single-collision J^Pt distribution Px_ col. applicable 
to W or Z production (Fig. 13"% 

Pi-co, (5» oc (5» a3 V2>Wi4,7. (38) 
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FIG. 38: The simulation (solid) and data (circles) distri- 
butions of (pi + u v ) for Z boson decays to muons (top) 
and electrons (bottom). The combined \ 2 is minimized 
in the fit for the recoil detector response parameters. 



The ^2 pt produced in a single minimum bias colli- 
sion can be different from that produced by the spec- 
tator partons in W or Z boson production. In order 
to allow for a difference, we scale the ^2pt drawn 
from the single-collision spectrum by a parameter 
Nw,z, which we tune on the Z boson data. 

With this model, the J2pt in a simulated event is 
obtained by adding the contributions from the spec- 
tator partons and the additional interactions. The 
corresponding recoil resolution is generated according 
to Eqn. I36l with a single tunable parameter N\y,z- 



3. Recoil Energy Resolution 

The measurement of the quark and gluon radiation 
is affected by detector energy resolution, which in 
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FIG. 39: The YltPr distribution for minimum bias data 
(solid) and a single pp collision (dashed), as derived from 
the minimum bias distribution. 



turn affects the measured recoil direction. We model 
the recoil angular resolution as a Gaussian distribu- 
tion with (70 = 0.14±0.01(stat), determined from fits 
to the A(/)(ut , —Px ) distribution in Z boson events 
(Fig. SOI)- Since the lepton directions are precisely 
measured, the width of the peak at A</> = is domi- 
nated by the recoil angular resolution. 

The energy resolution of the quark and gluon radi- 
ation is predominantly determined by stochastic fluc- 
tuations in the hadronic calorimeter, which motivate 
the functional form a UT oc \J u t rf ue . We measure the 
proportionality constant Shard using Z boson data. 

To tune Shard and Nw,z, we project the momen- 
tum imbalance Px + ux along the r\ and £ axes in 
Z boson decays (Fig. [44j) . The width of these pro- 
jections as a function of p l x provides information on 
N\y,z and Shard- At low the resolution is domi- 
nantly affected by Nw,z, with the Shard contribution 
increasing as the boson px increases. We compare 
the widths of the data and simulation projections as 
a function of p T and compute the x 2 . Minimizing 
this x 2 , we obtain N w ,z = 1.167 ± 0.026(stat) and 
Shard = [0.828 ±0.028(stat)] GeV 1 / 2 . The tuning is 
performed such that the statistical uncertainties on 
these parameters are uncorrelated. 



D. Recoil Model Cross-Checks 

The full recoil model, with parameters tuned from 
Z boson events, is applied to the simulated W boson 
sample. We compare the data to the predictions of 
distributions that can affect the final mass measure- 
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FIG. 40: The angle between the measured recoil and 
the direction opposite Px, for simulation (histogram) and 
data (circles) events where the Z boson decays to muons 
(top) or electrons (bottom). The x 2 from the Z — > fifi 
sample is minimized in the fit to the recoil angular resolu- 
tion. The corresponding uncertainty on raw is negligible. 



ment: the projections of the recoil along (u||) and 
perpendicular to (u±) the charged lepton; and the 
total recoil ut- 

The U|| distribution is directly affected by the mea- 
surements of lepton efficiency as a function of uii 
(Figs. HH and [TO)) and the modeling of lepton tower 
removal (Figs. [531 and |3"1| . The um is also sensitive 
to the boson px (Sec. HXBp and decay angular distri- 
butions, and to the recoil response and resolutions. 

Since ut is much less than the charged lepton px 
for our event selection, ft T ps \px + wJ. Thus, nix 
can be written as: 



rax ~ 2pTy 1 + u\\/pt ~ %Pt + u\\- (39) 
To a good approximation, any bias in u\\ directly en- 
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FIG. 42: The simulation (solid) and data (circles) um dis- 
tributions for W boson decays to \iv (top) and ev (bot- 
tom). The simulation uses parameters fit from Z boson 
data, and the uncertainty on the simulation is due to 
the statistical uncertainty on these parameters. The data 
mean (fj.) and RMS (a) are well-modeled by the simula- 
tion. 
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FIG. 41: The simulation (solid) and data (circles) p^ + u n 
(top) andp" + ?t{ (bottom) resolutions for Z boson decays 
to muons and electrons. The sum of the four \ 2 values is 
minimized in the fit for the recoil resolution parameters 
Nw,z an d Shard- Since there are four distributions and 
two fit parameters, each distribution contributes half a 
degree of freedom to the fit. 



ters as a bias in the wit fit. We compare the it|| 
distributions in data and simulation M|| in Fig. 1421 
and observe no evidence of a bias at the level of 
the data statistics and simulation systematics derived 
from the recoil model parameters. All backgrounds 
(Section lVIIip arc included in the comparison, except 
W — > tv, which has similar distributions to the other 
W leptonic decays. 

The u±_ distribution is dominantly affected by the 
recoil resolution, with a smaller contribution from the 
recoil response. The simulation models this distribu- 
tion well for both W — > ev and W — > [iv samples 
(Fig. H3). 

The mean of the ut distribution is sensitive to the 
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FIG. 43: The simulation (solid) and data (circles) u± dis- 
tributions for W boson decays to \w (top) and ev (bot- 
tom). The simulation uses parameters fit from Z boson 
data, and the uncertainty on the simulation is due to 
the statistical uncertainty on these parameters. The data 
mean (fj,) and RMS (a) are well-modeled by the simula- 
tion. 



FIG. 44: The simulation (solid) and data (circles) ut dis- 
tributions for W boson decays to \iv (top) and ev (bot- 
tom). The simulation uses parameters fit from Z boson 
data, and the uncertainty on the simulation is due to 
the statistical uncertainty on these parameters. The data 
mean (fi) and RMS (<r) are well-modeled by the simula- 
tion. 



recoil response and the boson p?, and is affected to 
a lesser extent by the resolution. The reverse is the 
case for the RMS of the ut distribution. Both are 
modeled well by the simulation for both W — > ev 
and W — > \iv samples (Fig. I31]) . 

The uncertainties on the mw fits from the recoil 
parameters (Table|V]) are determined by varying each 
parameter by ±3cr and assuming linear variation of 
the fit mw with the parameter. Since all uncertain- 
ties are uncorrelated, we add them in quadrature to 
obtain total recoil model uncertainties of 12, 17, and 
34 MeV on mw from the ttit, Pt, and ft T fits, re- 
spectively. The uncertainties are the same and 100% 
correlated for the electron and muon channels, since 



the recoil parameters are obtained from combined fits 
to Z — > ee and Z — > data. The uncertainty on 
the pt fit arises predominantly from the modeling of 
the ut < 15 GeV threshold used to select W boson 
events (Section IIV|) . 



VIII. BACKGROUNDS 

The event selection criteria (Section ITV|) result in 
W boson samples with high purity. However, the 
small residual backgrounds affect the distributions 
used for the mw fits. Both the W — > ev and W — > fiv 
samples receive contributions from: Z/-f* — > 11, 
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Input 



Shift (MeV) 



parameter 


rriT 


Pt 


Pt 


A 


-9 


-8 


2 


B 


-2 


15 


15 


N w ,z 


5 





22 


Shard 


5 


-3 


21 



TABLE V: Signed shifts in the mw fits due to la in- 
creases in the recoil model parameters. 



where one lepton is not detected; W — > tv, where the 
t decay products are reconstructed as a charged lep- 
ton; and multijet production, where at least one jet is 
misreconstructed. The W — > [iv sample also contains 
backgrounds from cosmic rays, where a muon passing 
through the COT is reconstructed on only one side of 
the COT, and long-lived hadrons decaying to /ivX, 
where the muon momentum is misreconstructed. 



A. W — > ev Backgrounds 
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FIG. 45: The f distribution of the W 



sample (triangles) and prediction (solid), with the f T 



ev candidate 
and 

m,T selection cuts removed. We fit for the normalization 
of the hadronic jet background (dotted) after fixing the 
normalization of the W — > ev distribution (dashed) in 
the peak. Not shown are the Z/7* — > ee and W — > rv 
backgrounds, whose relative normalizations are fixed from 
the simulation. 



We model the W — + tv and Z/^* — ► ee back- 
grounds using events generated with PYTHIA [5!| and 
simulated with a full GEANT-based detector simula- 
tion [13, SH . The full simulation models global de- 
tector inefficiencies and is thus more appropriate for 
predicting background normalizations than the cus- 
tom fast simulation. The multijet background is es- 
timated using a data-based approach. 

In the standard model the branching ratio for 
W — > ev is the same as for W — > tv, neglecting lep- 
ton masses. Measurements from LEP [19| test this 
prediction with a precision of 2.9%, and a slight dis- 
crepancy from the standard model is observed with 
a significance of 2.6o\ In estimating the W — » tv 
background, we assume the standard model predic- 
tion and determine the ratio of W — * tv to W — ► ev 
events from the ratio of acceptances of these two pro- 
cesses, as determined by the full GEANT-based detec- 
tor simulation. We include an uncertainty of 2.9%, 
corresponding to the statistical precision of the tests 
of this assumption. We estimate the W — > tv back- 
ground to be (0.93 ± 0.03)% of the W -> ev sample. 

The Z/j* background is determined from the ra- 
tio of Z/-f* — > ee to W — > ev acceptances deter- 
mined from the GEANT-based detector simulation, 
multiplied by the corresponding ratio of cross sec- 
tions times branching ratios. The ratio a ■ BR(Z — ► 
ee)/cr • BR(W — * ev) has been calculated in the stan- 
dard model to be 10.69 ± 0.08 |1 [H, and mea- 
surements are consistent with this value [26| [56| 



(65l |. We take an uncertainty of ±0.43 on this value 
from the CDF Run I measurement, and estimate the 
Z/7* — > ee background in the W — > ev candidate 
sample to be (0.24 ± 0.01)%. 

Multijet background enters the signal data sample 
when a hadronic jet is misreconstructed as an elec- 
tron and a second jet results in large ft T through en- 
ergy misreconstruction or the semi-leptonic decay of a 
hadron. To estimate this background, we remove the 
ft T threshold in our signal event selection to include 
the background-dominated kinematic region of low 
ft T . We then fit the observed f T spectrum to the com- 
bination of the hadronic jet, W — > ev, Z/j* — * ee, 
and W — > tv components, floating only the hadronic 
jet shape normalization (Fig. 145(1. 

In this fit, the shapes and normalizations for the 
W — > ev, Z/j* — * ee, and W — > tv components are 
determined from the GEANT-based simulation. The 
shape of the f T spectrum of the hadronic jet back- 
ground is determined from the single-electron events 
that pass an anti-electron identification requirement 
based on a neural network discriminant NN. The 
discriminant is determined by combining the electron 
quality variables (Section llVBp into a neural network 
[661 ] trained with single-electron data events, using 
f T to separate signal from background. 

Electron candidates in the W — > ev sample with 
low NN values have a high probability to be jets 
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misreconstructed as electrons. Events with such can- 
didates provide a p distribution characteristic of 
hadronic jet production. We apply a small correc- 
tion to this distribution to account for the expected 
contribution from W — > ev decay electrons with low 
NN values. 

This method relies on the assumption that the 
hadronic jet background has a p T distribution that is 
independent of the electron identification variables. 
As a test of this assumption, we perform the same fit 
for the jet background normalization, using only the 
isolation variable (Section llV A[) instead of the NN to 
select a hadronic jet subsample. We take a weighted 
average of the two fitted background normalizations, 
and assign an uncertainty that covers the range of 
the two results. The resulting background estimate 
is (0.25 ± 0.15)% of the W -> ev sample. 

The tot, Pt, an d p T distributions are obtained 
from the GEANT-based simulation for W and Z boson 
backgrounds, and from events in the W — > ev sam- 
ple with low- iViV electron candidates for the hadronic 
jet background. We fit these distributions (Fig. H5]l 
and include their shapes and relative normalizations 
in the mw template fits. The uncertainties on the 
background estimates result in uncertainties of 8, 9, 
and 7 MeV on mw from the tot, Pt, and p T fits, 
respectively (Table EJ). 



W-^ev backgrounds 





% of 




Srn w (MeV) 




Background 


W — > ev data 


m,T 


fit pr fit 


i T fit 


W ^ TV 


0.93 ± 0.03 


2 


2 


2 


Hadronic jets 


0.25 ± 0.15 


8 


9 


7 


Z/j* — > ee 


0.24 ± 0.01 


1 


1 





Total 


1.42 ± 0.15 


8 


9 


7 



TABLE VI: The percentages of the various backgrounds 
in the W — * ev data set, and the corresponding uncer- 
tainties on the m,T, pr and p fits for mw- 



B. W — + [w Backgrounds 

The W — > tv and Z/j* — > ///i backgrounds are 
modeled using events generated with PYTHIA [5!| and 
simulated with the GEANT [4(|-based detector simu- 
lation. We use the data to estimate backgrounds from 
cosmic rays, multijets, and hadrons decaying in flight 
to \xvX. 

Backgrounds from W — * tv and Z/'j* — ► fi[i to 
the W fiv sample are modeled in the same man- 
ner as for the W — > ev sample (Section I VIII A\ . We 
determine the ratio of the acceptance for W — > tv 
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FIG. 46: The parametrizations of the backgrounds to the 
W — > ev data sample. The backgrounds to the m,T (top), 
Pt (middle), and p (bottom) distributions are included 
in the mw fits. 
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or Z/-f* — > events to the acceptance for W — > /i^ 
events using the GEANT-based detector simulation. 
We assume equal branching ratios for the two W 
boson decay modes, and use the ratio a ■ BR(Z — > 
Wi)/a-BR{W -> H = 10.69±0.43 (Section|WIZ|)- 
We estimate the fraction of W tv (Z/j* — > /x/x) 
events in the W ^ [iv candidate sample to be 
(0.89 ±0.02)% [(6.6 ±0.3)%]. The Z/ 7 * -> p,p, back- 
ground is large because our event selection does not 
identify muons with \r\\ > 1.2. The tracker and muon 
detectors have incomplete or no coverage in the for- 
ward rapidity region, and the muons deposit only 
a few GeV of energy in the calorimeter. Thus, a 
Z/-f* — > fifj, event with one central and one forward 
muon is measured as a single-muon event with large 

Cosmic-ray muons passing close to the beam line 
are a source of background to the W \iv sample 
when the muon track is reconstructed on only one 
side of the COT. The cosmic-ray identification algo- 
rithm [57l | searches for unreconstructed tracks and re- 
moves cosmic rays with high efficiency. The residual 
cosmic-ray background is estimated using the recon- 
structed interaction time to and impact parameter 
do from the COT track fit. Figure [47l compares the 
to distributions of the W — > \w candidate sample, 
Z/j* — > fifi candidates, and identified cosmic rays. 
The cosmic ray fraction is fit by minimizing the % 2 of 
the sum of the Z/j* — > fxfj, and cosmic ray distribu- 
tions with respect to the W — ► fiv distribution. We 
obtain an alternative background estimate by com- 
paring the do distribution of identified cosmic rays to 
the do distribution of W — ► fj,v candidates with the 
do selection cut removed. The high impact parameter 
region of the W boson sample is enriched with cosmic 
rays, and is used to estimate the cosmic ray back- 
ground within the selection region |do| < 1 mm. We 
take the cosmic-ray background to be (0.05 ±0.05)%, 
where the uncertainty covers the range of results from 
the two estimates. 

Decay of a long-lived meson to a muon can re- 
sult in a reconstructed track with high momentum 
and large event f T . A low-momentum pion or kaon 
(< 10 GeV) that decays in the tracking chamber can 
be reconstructed as a high-momentum muon if the 
decay is in an azimuthal direction opposite the me- 
son's curvature (i.e., a kink in the trajectory). Such 
misreconstruction typically results in a poor COT 
track x 2 an d a large impact parameter. For each of 
these quantities we obtain a prompt muon distribu- 
tion from Z boson decays and a meson decay-in-flight 
distribution from the W boson sample by requiring 
either high COT track \ 2 or high impact parameter. 
We fit for the background fraction by summing the 




ip- i i U ill 

-15 -10 -5 5 10 15 

t (ns) 

FIG. 47: The to distributions for W — > candidates 
(solid), Z — > h/j, candidates (dashed), and identified cos- 
mic rays (dotted). For comparison, the distribution from 
Z candidates has been normalized to that of the W candi- 
dates, while the cosmic distribution has been normalized 
to 0.05% of the W boson sample. 



prompt muon distribution with the decay-in-flight 
distribution, and minimizing the x 2 with respect to 
the muon distribution from the W boson sample. We 
obtain a background fraction of (0.3 ± 0.2)%, where 
the uncertainty covers the range of the estimates ob- 
tained using the COT track \ 2 and impact parameter 
distributions. 

A separate class of hadronic background re- 
sults from high-momentum muons from short-lived 
hadronic decays, or energetic hadrons penetrating 
the calorimeter to the muon detectors. These back- 
ground muon candidates are typically accompanied 
by significant hadronic energy due to an associated 
hadronic jet, and can be separated using a muon iso- 
lation variable. Two such variables are determined 
by using either calorimeter energy or track momenta 
in an 77 — <fi cone of size 0.4 surrounding the muon 
candidate. Using the low f T region to select a jet- 
dominated sample, we fit the track and calorimeter 
isolation distributions of the W boson candidate sam- 
ple to the sum of the expected distributions from 
Z —>■ events and jet-dominated events. As a 
third method, we fit the ft T distribution, using muon 
candidates with high-isolation values to provide the 
f T distribution of the hadronic-jet background. From 
the range of results of the three methods, we obtain 
a jet background estimate of (0.1 ± 0.1)%. 

The distributions for the mw fit variables are ob- 
tained from the GEANT-based simulation for W and Z 
boson backgrounds, from identified cosmic ray events 
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for the cosmic ray background, and from events in 
the W —> [iv sample with high-% 2 (isolation) muons 
for the decay-in-flight (hadronic jet) background. In- 
cluding uncertainties on the shapes of the distribu- 
tions, the total uncertainties on the background es- 
timates result in uncertainties of 9, 19, and 11 MeV 



W->\iv backgrounds 



on mw from the mx, 
(Table En}- 



Pt, and p fits, respectively 





%of 




5m w (MeV) 




Background 


W — * fiv data 


TUT 


fit pt fit f 


> r fit 


z/i* -> m 


6.6 ± 0.3 


6 


11 


5 


W -^TV 


0.89 ± 0.02 


1 


7 


8 


Decays in flight 


0.3 ± 0.2 


5 


13 


3 


Hadronic jets 


0.1 ± 0.1 


2 


3 


4 


Cosmic rays 


0.05 ± 0.05 


2 


2 


1 


Total 


7.9 ± 0.4 


9 


19 


11 



TABLE VII: The percentages of the various backgrounds 
in the W — > \iv data set, and the corresponding uncer- 
tainties on the rriT, Pt and f T fits for mw- 



IX. PRODUCTION AND DECAY MODELS 

The measurement of the W boson mass relies on 
a complete model of W and Z boson production 
and decay. The production process is described by 
perturbative QCD and a parametrization of non- 
perturbative QCD effects, with parameters deter- 
mined from global fits to hadron-hadron and lepton- 
hadron collision data. W and Z boson decay are 
modeled using a next-to-leading-order electroweak 
calculation and includes QCD corrections for the lep- 
ton angular distributions, as a function of boson px- 
The most important process in the decay is photon 
radiation off the final-state charged lepton, which has 
been calculated at next-to- leading order [621 ]. 



A. Parton Distribution Functions 

The longitudinal momentum of the produced W 
or Z boson depends on the momenta of the inter- 
acting partons. These momenta, generally expressed 
in terms of the fractions x\ of the colliding (anti- 
proton energies, are not known on an event-by-event 
basis. The Xi parton distribution functions (PDFs) 
are however well constrained by hadron-hadron and 
lepton-hadron collision data. The distributions have 
been parametrized as simple functional forms for the 
quarks, antiquarks, and gluons inside a proton. Two 
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FIG. 48: The parametrizations of the backgrounds to the 
W — > \iv data sample. The backgrounds to the tut (top), 
Pt (middle), and p (bottom) distributions are included 
in the mw fits. Not shown are the small hadronic-jet and 
cosmic-ray background distributions. 
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independent fits to the glob al data, performed by the 
MRST [23] and CTEQ [HI collaborations, constrain 
the parameters in these PDFs. 

We model the quark momentum fractions using the 
next-to-leading-order CTEQ6M parton distribution 
functions. The CTEQ parametrization [24| for most 
of the distribution functions inside the proton is: 

x P fa(x p , Q ) = A x At (1 - x p ) A2 e A3X "(l + A iXp ) M , 

(40) 

where f a are the distributions of a particular quark or 
gluon combination a, A% are the fit parameters, and 
Qo is the energy scale at which the parameters are 
defined. The functions at a particular energy scale Q 
are determined by a perturbative evolution calcula- 
tion known as the DGLAP equation 67]. 

The uncertainty on the mw measurement arising 
from uncertainties on the PDF parameters is de- 
termined using a set of 40 PDFs provided by the 
CTEQ collaboration. The set covers the ±1.6<r 
(90% C.L.) uncertainties (68| for the eigenvectors of 
the parametrization. The mass shift of a particu- 
lar + 1.6er PDF, relative to the corresponding —1.6<r 
PDF, determines the uncertainty due to that eigen- 
vector. We calculate the total PDF uncertainty using 
the quadrature sum of all eigenvector contributions 
El: 



5mPpF = re 



\rrf+-m i w ) 2 



(41) 



where m 1 ^ represents the mass fits for the ±1.6(7 
shifts in eigenvector i. These fits are performed us- 
ing templates and simulated pseudoexperiments both 
generated with pythia [59]. The resulting 5m w DF 
are 11, 9, and 13 MeV, for the mx, Pt, and f T fits, 
respectively. A fit to pseudodata using the MRST 
PDF set results in mw shifts smaller than these un- 
certainties. 



B. W and Z Boson pr 

Because mass is a Lorentz invariant, the W boson 
transverse mass is only weakly sensitive to the W 
boson transverse momentum p^ . However, the decay 
lepton pt spectra are more significantly affected by 
the pjf distribution. 

At the Tevatron, the pt spectra of W and Z bosons 
peak at a few GeV (Fig. 25]), where the shapes are 
predominantly determined by non-perturbative QCD 
interactions. We model the distribution with the 
RESBOS generator [6111 . which uses the Collins-Soper- 
Sterman (CSS) [69| resummation formalism and a 



parametrized non-perturbative form factor. In this 
formalism, the cross section for W boson production 
is written as: 



da(pp W + X) _ 1 
2^ 2 



dsd 2 prp dy 



= 7T5${s-m 2 w ) 



x W fk (b,s,x % ) + Y{p T J ,s,Xi), 



(42) 



where the parton energy fractions of the (anti- 

proton, y = 0.5 hi{xp/xp) is the boson rapidity, b is 
the relative impact parameter of the partons in the 
collision, Y is a function calculable at fixed order, 
and W can be separated into its perturbative and 
non-perturbative components. We use the Brock- 
Landry-Nadolsky-Yuan (BLNY) form for the non- 
perturbative component: 



e [si-32 ln(5^)-gifl> 3 M^00x p x p )]b 2 



where Q = 1.6 GeV and gt are parameters suggested 
by the CSS formalism to be universal to processes 
with initial state quarks and colorless objects in the 
final state [oil ]. 

The <?2 parameter affects the position of the most 
probable p^ and is the most relevant for the mw 
measurement. We use g\ = 0.21 GeV 2 , gi = 0.68 
GeV 2 , and g% = —0.60, which are determined from 
fits to global Drell-Yan data [6l|. We verify the ap- 
plicability of these values to our data by fitting the 
dilepton px distribution (Fig. |4T))) for g2- We find 
g 2 = [0.685 ± O.Q48(stat)] GeV 2 , consistent with the 
global fits. Varying gi by ±3cr in pseudoexperiments 
and taking the fit raw to be linearly dependent on 
g2, we find that the uncertainty of 5g 2 = 0.048 GeV 2 
results in uncertainties of 3, 9, and 5 MeV, on mw for 
the mr, Pt, and f T fits, respectively. These uncer- 
tainties are the same and 100% correlated between 
the electron and muon channels, since gi is fit us- 
ing Z — > ee and Z — * /i/i. Neglecting correlations 
between PDFs and the gi parameters, we find that 
uncertainties on the other gi do not significantly af- 
fect the mw measurement. 



C. W Boson Decay 

The mw measurement is sensitive to the charged 
lepton decay angle relative to the boson px- The 
mismodcling of this angle can bias the projection of 
the recoil along the lepton (un), which in turn affects 
mw measured from the mj fit (Section IVIIDp . 

The lepton decay angle is predicted by the matrix 
element calculation in the RESBOS generator, which 
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tion can be written as [711 ] . 
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FIG. 49: The simulation (solid) and data (circles) 
distributions for Z boson decays to fifi (top) and ee (bot- 
tom). The distributions are used to fit for the non- 
pcrturbative parameter Q2, which determines the most 
probable value of pf. Since there are two distributions 
for one fit parameter, each distribution contributes half 
of a degree of freedom. The mean (p) and RMS (a) are 
consistent between data and simulation. 



computes the differential cross section dp w d y^ d Q2 dn - 
The angular distributions are defined in the Collins- 
Soper rest frame of the W boson [70| . In this frame, 
the z-axis is defined to bisect the angle between the 
proton momentum and the opposite of the antiproton 
momentum. 



The angular component of the differential cross sec- 



^ oc(l + cos 2 6) + i A (l - 3 cos 2 6) 

+ Ai sin 29 cos 4> + - A 2 sin 2 9 cos 2(f> 

+ A3 sin 9 cos <f> + A± cos 9 + A$ sin 2 9 sin 2<p 
+ A 6 sin 20 sin (/> + Aysinflsin^, (44) 

where the Ai(p}p,yw) have been calculated to NLO 
in a s . Because of the V — A structure of the elec- 
troweak interaction, for leading-order valence quark 
interactions all Ai are zero except A4 = 2. The 
Ai can be determined experimentally through a mo- 
ments analysis [72| of the lepton angle in the Collins- 
Soper frame. 

We have performed a moments analysis to extract 
the Ai from the resbos generator, and compared the 
results to those obtained [72j from the DYRAD event 
generator [73j |. which produces W+ jet events to or- 
der a 2 . The two generators give consistent results in 
the overlapping region 15 GeV < < 100 GeV. 



D. Photon Radiation 

The quarks, the W boson, and charged lepton have 
non-zero electromagnetic charge and can radiate pho- 
tons in the W boson production process. Radiation 
off the initial-state quarks and the W boson propa- 
gator have a negligible effect on the invariant mass 
distribution of the W boson. Radiation off the final- 
state charged lepton reduces the measured transverse 
mass (relative to the W boson mass) and must be ac- 
curately modelled. 

We study photon radiation using the WGRAD event 
generator [62| , which models the full next-to- leading- 
order (NLO) electroweak physics. The generator al- 
lows an independent study of photon radiation from 
the initial-state quarks (ISR) , the W boson propaga- 
tor, and the final-state charged lepton (FSR). Inter- 
ference between the contributing diagrams can also 
be studied independently. We verify that the initial- 
state, propagator, and interference effects do not af- 
fect the measured W boson mass, within the 5 MeV 
statistical uncertainty of the simulation. 

We simulate final-state photon radiation in our 
RESBOS-generated W and Z boson events by gener- 
ating a photon for each charged lepton. The energy 
and angular spectra are taken from the WGRAD gen- 
erator using the appropriate boson mass. To avoid 
the infrared divergence that arises when the photon 
momentum goes to zero, we require E 1 > 5 s ys/2, 
where S s = 10~ 4 . We find that increasing 5 S to 10 -3 
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FIG. 50: The cube root of the fraction of electron (solid) 
or muon (dashed) momentum contained in the radiated 
photon (top), and the square root of the angle A_R (bot- 
tom) between the radiated photon and the electron (solid) 
or muon (dashed). 



does not affect the mw measurement, at the level of 
the 5 MeV statistical precision of the tests. 

The energy of a photon in a given event is cal- 
culated from the fraction y = E^/Ei. The photon 
angle AR = \J (Ar;) 2 + (Acf>) 2 is taken with respect 
to the charged lepton. To improve the phase space 
sampling, we sample from a two-dimensional distri- 
bution of the variables V AR and y 1 / 3 when selecting 
a photon. The individual distributions of these vari- 
ables are shown in Fig. l50l 

We validate our photon simulation by fitting a sam- 
ple of events generated with FSR using WGRAD to 
templates generated with leading-order WGRAD and 
photons simulated according to our model. We find 
our FSR model to be consistent with that of WGRAD 
at the level of the 5 MeV statistical precision of the 



Fit 


5m w ,z(n) (MeV) 


5mw,z(e) (MeV) 


UlT 


-158 


-138 


PT 


-206 


-186 


f T 


-77 


-59 


m n 


-196 


-215 



TABLE VIII: The mass shifts obtained by fitting events 
generated with our simulation of single-photon radiation 
to templates generated without final-state photon radia- 
tion. The shifts are for the W boson mr, Pt, and ft T fits, 
and for the Z boson ran fit. The shifts have statistical 
uncertainties of 7 MeV each. 



test. The total effect of including FSR is shown in 
Table IVIIII Since FSR reduces the charged lepton 
momentum, the shift is largest for the px fit and 
smallest for the f T fit. The effects are smaller for elec- 
trons than for muons because the electron calorimeter 
energy measurement recovers much of the energy of 
FSR photons. 

We approximate the effect of next-to-next-to- 
leading-order FSR by increasing the photon's mo- 
mentum fraction (y) by 10%, consistent with the re- 
sults of a study of higher-order photon radiation [74( . 
We take half the correction as a systematic uncer- 
tainty to account for higher-order QED effects. 

The total uncertainty due to photon radiation is 
the quadrature sum of: uncertainties on ISR, inter- 
ference between ISR and FSR, and radiation off the 
propagator (5 MeV); uncertainty due to the infrared 
cutoff of the FSR photon (5 MeV); the FSR model 
(5 MeV) ; and uncertainties on higher-order FSR cor- 
rections (7 MeV for the electron and 8 MeV for the 
muon m,T fits). The total uncertainties are 12 (11), 
13 (13), and 10 (9) MeV, for the muon (electron) tot, 
Pt, and f T fits, respectively. 



X. W BOSON MASS FITS 

We fit the W boson data distributions to a sum 
of background and simulated signal templates of the 
wit, Pt, and ^distributions, fixing the normaliza- 
tion of the sum to the number of data events. The fit 
minimizes the negative log likelihood (Section III Dp 
as a function of the template parameter mw, which 
is defined by the relativistic Breit-Wigner mass dis- 
tribution [ill ]: 



da 
dm 



(to 2 - to 2 ^) 2 + toT 2 ^/ 



2 ' 



(45) 



where m is the invariant mass of the propagator. The 
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likelihood is calculated in mw steps of 1 MeV. We use 
the standard model W boson width Tyy = 2.094 GeV, 
which has an accuracy of 2 MeV and is calculated 
for mw = 80.393 GeV. Using pseudoexperiments, we 
find the input T\y affects the fit mw according to the 
relation dmw/dTw = 0.14 ± 0.04. 



A. Fit Results 

The results of the mr fits are shown in Fig. I5T1 
and Table [LXl gives a summary of the 68% confidence 
level uncertainties associated with the fits. We fit for 
mw in the range 65 GeV < < 90 GeV, where 
the fit range has been chosen to minimize the total 
uncertainty on mw- The pt and f T distributions are 
fit in the range 32 GeV < p T < 48 GeV (Figs. 
and 1531 respectively) and have uncertainties shown 
in Tables [X] and IXI1 respectively. We show the in- 
dividual fit results in Table IXII1 and the negative 
log-likelihoods of all fits in Fig. [54l 

tut Fit Uncertainties 
Source W — * fiv W — ► ev Correlation 



Tracker Momentum Scale 


17 


17 


100% 


Calorimeter Energy Scale 





25 


0% 


Lepton Resolution 


3 


9 


0% 


Lepton Efficiency 


1 


3 


0% 


Lepton Tower Removal 


5 


8 


100% 


Recoil Scale 


9 


9 


100% 


Recoil Resolution 


7 


7 


100% 


Backgrounds 


9 


8 


0% 


PDFs 


11 


11 


100% 


W Boson pt 


3 


3 


100% 


Photon Radiation 


12 


11 


100% 


Statistical 


54 


48 


0% 


Total 


60 


62 




TABLE IX: Uncertainties 
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FIG. 51: The simulation (solid) and data (points) mx dis- 
tributions for W boson decays to pv (top) and ev (bot- 
tom). The simulation corresponds to the best-fit mw, 
determined using events between the two arrows. The un- 
certainty is statistical only. The large x 2 f° r the electron 
fit is due to individual bin fluctuations (Fig. I55|l and does 
not bias the fit result, as evidenced by the small change 
in the fit mw when the fit window is varied (Fig. 158 p . 



to combine all six mw fits. When combining any sub- 
set of fits, the appropriate smaller covariance matrix 
is used. 

The result of combining the mw fits to the mj 
distribution in the W — > fiv and W — * ev channels is 



We combine results from the W — + fiv and W — * 
ev fits u sing the Best Linear Unbiased Estimator 
(BLUE) [tH| . The BLUE algorithm defines a pro- 
cedure for constructing a complete covariance ma- 
trix using the derivative of mw with respect to each 
model parameter [l8j]. We construct this matrix as- 
suming each source of systematic uncertainty is in- 
dependent of any other source of uncertainty. The 
resulting covariance matrix (Table IXIIIj) is then used 



m w = 80.417 ± 0.048 GeV. (46) 

The x 2 /dof of the combination is 3.2/1 and the prob- 
ability that two measurements of the same quantity 
would have a x 2 /dof at least as large as this is 7%. 

The combination of the fits to the px distribution 
yields 

m w = 80.388 ± 0.059 GeV, (47) 
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FIG. 52: The simulation (solid) and data (points) pr dis- 
tributions for W boson decays to \iv (top) and ev (bot- 
tom). The simulation corresponds to the best-fit mw, 
determined using events between the two arrows. The 
uncertainty is statistical only. 



with a x 2 /dof of 1.8/1 and an 18% probability for the 
two measurements to obtain a % 2 /dof > 1.8. 

The results of the fits to the f T distribution gives 

m w = 80.434 ± 0.065 GeV, (48) 

with a 43% probability of obtaining a x 2 /dof at least 
as large as observed (0.6/1). 

Combining the rriT, Pt, and f T fits 
within the individual decay channels gives 
m w = (80.352 ±0.060) GeV with a x 2 /dof of 1.4/2 
for the W -> [iv channel and m w = (80.477 ± 0.062) 
GeV with a x 2 /dof of 0.8/2 for the W — > ev channel. 

We combine the six fits with the BLUE procedure 
to obtain our final result of 



m w = 80.413 ± 0.048 GeV, 



(49) 
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TABLE X: Uncertainties in units of MeV on the charged 


lepton transverse momentum fit for mw in 


the W — > (j,v 


and W — ► ev samples. 








f T Fit Uncertainties 
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which has statistical and systematic uncertainties of 
34 MeV each. The statistical correlations between 



TABLE XI: Uncertainties in units of MeV on the missing 
transverse momentum fit for mw in the W ~ > [iv and 
W — * ev samples. 



the fits, determined from simulation pseudoexperi- 
ments, are shown in Table IXIVl The relative weights 
of the fits are 47.7% (32.3%), 3.4% (8.9%), 0.9% 
(6.8%) for the ray, Pt and ■p fit distributions, re- 
spectively, in the muon (electron) channel. The com- 
bination establishes an a priori procedure to incor- 
porate all the information from individual fits, and 
yields a x 2 /dof of 4.8/5. The probability to obtain a 
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FIG. 54: The negative log of the likelihood ratio C/Co, 
where Co is the maximum likelihood, as a function of mw 
for the rriT (solid), pr (dashed), and p T (dotted) fits in 
the electron (top) and muon (bottom) channels. 



FIG. 53: The simulation (solid) and data (points) 
p T distributions for W boson decays to fiu (top) and 
ev (bottom). The simulation corresponds to the best- 
fit mw, determined using events between the two arrows. 
The uncertainty is statistical only. 



Distribution 


mw (GeV) 


X 2 /dof 


mr(e, v) 


80.493 ± 0.048 ± 0.039 


86/48 


p T (e) 


80.451 ±0.058 ±0.045 


63/62 


P T (e) 


80.473 ± 0.057 ± 0.054 


63/62 


mrin, v) 


80.349 ± 0.054 ± 0.027 


59/48 




80.321 ±0.066 ±0.040 


72/62 


i> T (a*) 


80.396 ± 0.066 ± 0.046 


44/62 



TABLE XII: The results of the fits for mw to the mr, 
PT, and p distributions in the electron and muon decay 
channels. The first uncertainty is statistical and the sec- 
ond is systematic. 



m T (e,v) m T (p>,,v) Pr(e) Pt{^) jL (g) p T (p) 



mT{e,u) 64 2 

mr(lA, v) 

p T (e) 

pr(ji) 

t T (e) 

* T (M) 



27 2 
61 2 



61 2 


27 2 


61 2 


28 2 


27 2 


59 2 


28 2 


59 2 


75 2 


35 2 


51 2 


32 2 




77 2 


32 2 


53 2 






81 2 


43 2 








81 2 



TABLE XIII: The complete covariance matrix for the mr, 
Pt, and p T fits in the electron and muon decay channels, 
in units of MeV 2 . The matrix is symmetric. 



X 2 at least as large as this is 44%. 



B. Cross-Checks 



Figures l55ll57l show the differences between data 
and simulation, divided by the statistical uncer- 
tainties on the predictions, for the mr, Pt and 
p T distributions. Figures I58ll60l show the variations 



49 



Correlation 


W -> fiu 


W -> ev 


rriT — Pt 


70% 


68% 


rrvr—f T 


72% 


63% 


VT~f T 


38% 


17% 



TABLE XIV: The statistical correlations between the 
rriT, pr, and p fits in the electron and muon decay chan- 
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FIG. 55: The difference between the data and simulation, 
divided by the statistical uncertainty on the prediction, 
for the tut distributions in the muon (top) and electron 
(bottom) channels. The arrows indicate the fit region. 
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FIG. 56: The difference between the data and simulation, 
divided by the statistical uncertainty on the prediction, 
for the pr distributions in the muon (top) and electron 
(bottom) channels. The arrows indicate the fit region. 



The largest systematic uncertainty in this difference 
arises in the muon channel from the uncertainty on 
the alignment parameters ao and a 2 (Table [TTT|) . The 
uncertainties on the mass difference due to these pa- 
rameters are 49 MeV and 56 MeV, respectively, for 
a total uncertainty of 75 MeV. Any bias in these pa- 
rameters affects the positive and negative lepton fits 
in opposite directions, and thus has a negligible net 
effect when the two are combined. 



of the fitted mass values, relative to the nominal re- 
sults, as the fit regions are varied. These plots show 
variations consistent with statistical fluctuations. 

The variation of the px fits with time, detector re- 
gion, and lepton charge ( Table |XV|) show no evidence 
of dependence on time or detector region. There 
is a difference between positive and negative lepton 
mass fits at the level of » 2er in each decay channel. 



XI. SUMMARY 

We have performed a measurement of the W boson 
mass using 200 pb _1 of data collected by the CDF II 
detector at y^s = 1.96 TeV. From fits to px, and 
fi T distributions of the W — ► fiv and W — > ev data 
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FIG. 57: The difference between the data and simulation, 
divided by the statistical uncertainty on the data points, 
for the ft T distributions in the muon (top) and electron 
(bottom) channels. The arrows indicate the fit region. 



Fit difference W — > \w W — > ev 

m w {l + )-mw{l~) 286 ± 152 257 ± 117 

m w (<j)i > 0)-m w (<f>i < 0) ± 133 116 ± 117 
mw(Mar, 2002- Apr, 2003) - 

m^(Apr, 2003-Sep, 2003) 75 ± 135 -107 ± 117 



TABLE XV: Differences of mw in the pr fits between 
positively and negatively charged leptons, leptons in the 
upper and lower halves of the detector, and early and late 
data. The units are MeV. 



samples, we obtain 

mw = 80.413 ± 0.048 GeV, (50) 

which is the single most precise determination of mw 
to date. The uncertainty includes statistical and sys- 
tematic contributions of 34 MeV each. 



FIG. 58: The shifts in measured mw for variations in 
the lower (top) and upper (bottom) edges of the wit fit 
range. The electron channel is denoted by open triangles 
and the muon channel by solid circles. The error bars in- 
dicate the expected statistical variations from simulation 
pseudoexperiments. The dashed (dotted) lines indicate 
the statistical uncertainty on the mw fit using the de- 
fault fit range in the electron (muon) channel. 



Combining this result with the Run I Tevatron 
measurements using the method in [76| . we obtain 
a CDF Run I/II combined result of 

m w = 80.418 ± 0.042 GeV, (51) 

and a combined Tevatron result of 

m w = 80.429 ± 0.039 GeV. (52) 

In these combinations, we take the uncertainties due 
to PDFs and photon radiation to be fully correlated 
between our measurement and the previous Tevatron 
measurements. In the BLUE combination method 
[75| . each uncertainty source contributes its covari- 
ance matrix, and all covariance matrices are summed 
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FIG. 59: The shifts in measured mw for variations in 
the lower (top) and upper (bottom) edges of the pr fit 
range. The electron channel is denoted by open triangles 
and the muon channel by solid circles. The error bars in- 
dicate the expected statistical variations from simulation 
pseudoexperiments. The dashed (dotted) lines indicate 
the statistical uncertainty on the mw fit using the de- 
fault fit range in the electron (muon) channel. 



FIG. 60: The shifts in measured mw for variations in 
the lower (top) and upper (bottom) edges of the f T fit 
range. The electron channel is denoted by open triangles 
and the muon channel by solid circles. The error bars in- 
dicate the expected statistical variations from simulation 
pseudoexperiments. The dashed (dotted) lines indicate 
the statistical uncertainty on the mw fit using the de- 
fault fit range in the electron (muon) channel. 



to obtain the total covariance matrix. We evaluate 
an individual contribution to the uncertainty on our 
result by ignoring its respective covariance matrix 
and repeating the six-fold combination of our individ- 
ual electron and muon channel tot, Pt and f T fits. 
The difference in quadrature between the total un- 
certainty, including and excluding a given covariance 
matrix contribution, is taken to be the uncertainty 
due to that source. Following this procedure, we ob- 
tain the systematic uncertainty contributions due to 
PDFs and QED radiative corrections to be 12.6 MeV 
and 11.6 MeV respectively (Table IXVip . for a com- 
bined uncertainty of 17.2 MeV. 

Assuming no correlation between the Tevatron and 



LEP measurements, we obtain a new world average 
of 

m w = 80.398 ± 0.025 GeV. (53) 

Our measurement reduces the world uncertainty to 
31 parts in 10 5 , and further constrains the properties 
of the Higgs boson and other new particles coupling 
to the W and Z bosons. Within the context of the 
standard model, fits made to high energy precision 
electroweak data in 2006 gave tuh = 85+28 GeV, with 
m H < 166 GeV at the 95% confidence level [Hj]. The 
values used for the top quark and W boson masses 
in these fits were m t = (171.4 ± 2.1) GeV and mw = 
(80.392 ± 0.029) GeV, respectively. Updating these 
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Source 



Uncertainty (MeV) 



Lepton Scale 23.1 

Lepton Resolution 4.4 

Lepton Efficiency 1.7 

Lepton Tower Removal 6.3 

Recoil Energy Scale 8.3 

Recoil Energy Resolution 9.6 

Backgrounds 6.4 

PDFs 12.6 

W Boson p T 3.9 

Photon Radiation 11.6 

TABLE XVI: Systematic uncertainties in units of MeV on 
the combination of the six fits in the electron and muon 
channels. Each uncertainty has been estimated by remov- 
ing its covariance and repeating the six- fold combination. 



fits with the most recent world average values of mt = 
(170.9 ± 1.8) GeV and m w = (80.398 ± 0.025) GeV 
[Eg. (1531) | , and using the methods and data described 
in [Tg| and $7% , gives m H = 76+24 GeV ; with m H < 



144 GeV at the 95% confidence level. The effect of 
the new mw value alone is to reduce the predicted 
value of the standard model Higgs boson mass by 6 
GeV. 

We anticipate a significant reduction in the uncer- 
tainty of future CDF mw measurements using larger 
available data sets. The dominant uncertainties on 
this measurement are due to W boson statistics and 
to the lepton energy scale calibration (Table IXV1[) , 
and will be reduced with increased statistics in the 
W boson and calibration data samples. 



APPENDIX A: ELECTRON AND PHOTON 
INTERACTIONS 



The simulation of electrons and photons (Sec- 
tion If 11 B 2[) uses the Bethe-Heitler differential cross 
sections for electron bremsstrahlung and photon con- 
version [48|. Defining y as the final state energy di- 
vided by the initial state energy, the bremsstrahlung 
cross section is: 



da 
dy 



ka.EMr e 



4 

3y 



y )MZ) 



r 



(Al) 



where: 

ipi(Z) = Z 2 [ln(184.15Z- 1 /3) _ fl + 

Zln(1194Z- 2 / 3 ), 
MZ) = Z 2 + Z, (A2) 
/ = a 2 [(l + a 2 )- 1 +0.20206 - 0.0369a 2 + 
0.0083a 4 - 0.002a 6 ], 

and a = olemZ. We define a material's radiation 
length Xo according to jUj : 



An 



4a E Mr 2 e N A (npi(Z)/A, 



(A3) 



where p is the density of the material. In terms of 
the radiation length, the cross section is 



da 



A 



dy N A X p 
where 



y 



G 



MZ)/9 

MZ) ' 



(A4) 



(A5) 



The conversion cross section takes a similar form, 
since the relevant Feynman diagram is a rotation of 
the bremsstrahlung process [H| : 

T- = ^7 J Y- t 1 ^ ( 4 /3 + C)y(l - y)\ . (A6) 
dy N A X p 

The Compton scattering cross section as a func- 
tion of scattering angle is given by the Klein- Nishina 
formula (H: 



da 
dtt 



n 

2 



1 



cos 



k 2 (l 



COS ( 



[1 + fc(l - cos 0)] 2 [1 + jfe(l - COS t 



w. 

(A7) 

where k = E y /m e . The scattering angle is kinemati- 
cally related to the energy loss by [54[ 

y = k'/k= [l + /s(l- cos 6>)n\ (A8) 

where k' is the energy of the photon after scattering, 
in units of m e . Using this equation, the differential 
cross section with respect to y can be written as [4fj| | : 

da ... 
— cx 1 y + y, A9 
dy 
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ignoring terms containing 1/k. 
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